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ABSTRACT 

Problems  of  combustion  efficiency,  incomplete  prevaporization 
and  flame  stabilization  (blowoff  and/or  flashback)  associated  with 
advanced  prevaporizing/premixing  turbojet  combustors,  afterburners 
and  ramjet  combustors  are  examined  experimentally  using  a sinplified 
axisymmetric  burner.  This  flameholder  configuration  retains  the  same 
fundamental  combustion  processes  of  prevapori zing/premixing  combustors 
while  eliminating  some  of  the  complexity  of  the  practical  hardware. 

In  addition  to  variations  in  combustor  pressure,  inlet  tenperature, 
equivalence  ratio  and  air  velocity  typical  of  modem  combustors, 
geometry  and  heterogeneous  effects  are  also  considered. 

Both  flame  stabilization  and  conbustion  efficiency  data 
are  examined  using  conventional  (Lefebvre,  1966)  and  characteristic 
time  (Mellor,  1976)  correlations.  Characteristic  times,  which  have 
already  been  used  to  correlate  and  predict  gaseous  emissions  in 
conventional  conbustors,  are  quantified  here  for  blowoff.  The  model 
is  based  on  the  flame  stabilization  theory  of  Zukoski  and  Marble  (1956) 
and  includes  inlet,  geometry  and  heterogeneous  effects  simultaneously, 
unlike  traditional  loading  parameters.  Conbustion  efficiency  data 
obtained  for  this  study  are  not  yet  included  in  the  correlation. 
Detailed  internal  species  concentration,  tenperature  and  velocity 
measurements  are  planned  to  help  identify  proper  scaling  parameters 
for  the  conbustion  efficiency  model. 

Flame  flashback  into  the  fuel  preparation  tube  has  not  been 
observed  for  any  of  our  operating  conditions.  These  results  are 
consistent  with  a recent  literature  review  (Plee  and  Mellor,  1977) 
which  has  identified  important  mechanisms  of  tpstream  flame  propaga- 
tion. This  study  shows  that  flow  disturbances  rather  than  conbustor 
inlet  conditions  are  responsible  for  the  "flashback"  reported  in 
non- catalytic  conbustors. 
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I . INTRODUCriCW  AND  SUNWARY 


Gas  turbine  combiostion  systems  currently  provide  the  major 
source  of  propulsion  for  both  military  and  commercial  aircraft. 

In  addition,  continuous  flow  combustors  are  used  in  ramjets  and  are 
becoming  increasingly  iji^x)rtant  as  a secondary  energy  source  in 
stationary  power  plant  applications.  In  the  latter  case,  they  provide 
an  effective  method  for  satisfying  peak  load  requirements.  These 
combustors,  many  of  \diich  were  originally  developed  for  aircraft 
applications,  are  designed  to  maximize  the  amount  of  energy  release  per 
unit  volume.  Fig.  1-1  illustrates  schematically  a typical  gas  turbine 
conijustor  consisting  of  three  zones,  primary  combustion,  secondary 
combustion  and  dilution,  which  are  usually  established  by  discrete 
air  addition  through  the  liner  wall. 

The  primary  zone  normally  extends  from  the  injector  to  the 
first  set  of  air  addition  holes.  This  region  is  characterized  by  fuel 
vaporization  and  recirculating  air  flow  as  well  as  the  most  intense 
combustion.  Flame  stabilization  is  acconplished  by  swirling  the  inlet 
flow,  air  penetration  jet  inpingement  and  sudden  expansion  of  the 
combustor  liner.  The  creation  of  hot  recirculating  flow  regions 
provides  heat  and  free  radicals  to  the  shear  layer  and  serves  as  a 
continuous  ignition  source;  the  shear  layer  is  a region  of  extremely 
high  turbulence  located  between  the  low  velocity  central  recirculation 
zone  and  the  high  wlocity  inlet  air  stream.  Liquid  fuel  injected  into 
the  shear  layer  rapidly  vaporizes  and  reacts  with  the  incoming  air 
forming  a hollow  shear  layer  flame  zone  which  surrounds  the  central 
recirculation  zone  and  eventually  closes  to  the  centerline  (Mellor 
et  al.,  1972;  Tuttle  et  al.,  1973).  Tuttle  et  al.  (1973)  and  Shisler 
et  al.  (1974)  have  probed  these  flames  and  found  that  the  hollow  shear 
layer  flame  resembles  a turbulent  diffusion  flame.  Because  of  the  very 
nature  of  this  type  of  flame  and  the  fact  that  very  little  air  is 
added  in  the  primary  zone,  both  local  and  overall  equivalence  ratios 
are  high  providing  relatively  wide  ignition  and  blowoff  limits. 
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The  combustion  process  continues  to  completion  in  the 
secondary  zone.  Here  additional  air  is  provided  to  further  oxidize 
the  remaining  carbon  monoxide  and  unbumed  hydrocarbons.  Finally,  the 
products  of  combustion  are  further  cooled  in  the  dilution  zone.  It  is 
the  purpose  of  this  section  to  lovfer  the  exhaust  gases  to  an  acceptable 
tenperature  for  the  turbine  inlet.  Thus,  the  overall  combustor 
operates  quite  fuel  lean;  however,  locally  high  equivalence  ratios 
and  temperatures  are  present  in  the  primary  zone.  Typical  combustion 
efficiencies  for  this  type  of  burner  are  in  excess  of  99%  at  full  power. 

Although  local  stoichiometric  regions  provide  excellent 
ignition  and  blowoff  characteristics,  they  are  also  the  primary 
source  of  NO^  since  extremely  high  temperatures  are  present.  Mellor 
(1976)  and  Verkamp  et  al . (1974)  have  compared  the  emission  levels 
from  conventional  combustors  with  the  proposed  aircraft  emission  goals 
of  the  future  (EPA,  1973;  Blazowski  and  Henderson,  1974)  and  concluded 
that  major  burner  modifications  are  required  to  meet  these  new  stan- 
dards. The  problem  is  illustrated  in  Fig.  1-2  for  the  EPA  class  T-2 
regulations.  Conventional  conbustors  do  not  appear  capable  of  meeting 
both  the  NO^  and  CO  emission  levels  simultaneously  over  their  entire 
operating  range;  instead,  CO  quenched  by  the  cold  incoming  air  and 
NO^  formed  in  high  tenperature  regions  of  the  flow  are  problems  at  idle 
and  maximum  power,  respectively.  Minor  burner  modifications  to  con- 
ventional combustors  only  allow  movement  within  the  shaded  region 
rather  than  in  a direction  toward  the  origin  (Mellor,  1976;  Verkamp 
et  al.,  1974). 

In  an  effort  to  simultaneously  reduce  NO^  and  CO  levels, 
advanced  combustor  concepts  involving  fuel  prevaporization  and  premixing 
are  currently  being  considered  for  aircraft  (Roffe  and  Ferri,  1975  and 
1976;  Roberts,  P.  et  al.,  1976;  Niedzwiecki,  1975;  Fear,  1976;  Roberts, 
R.  et  al.,  1975  and  1976;  Mosier  and  Roberts,  1974;  Bahr,  1975;  Gleason 
et  al.,  1976;  Roffe,  1976),  automotive  (Spadaccini  and  Szetela,  1975; 
Azelbom  et  al.,  1973),  ramjet  (Stull  et  al.,  1974)  and  catalytic 
(Anderson  et  al.,  1975;  Blazowski  and  Bresowar,  1974;  Blazowski  and 
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Conventional 

Combustors 


1979  EPA  Class 
T2  Regulations 
SFC  = 0.35 
(lb  fuel/ lb 
thrust  hr) 


Advanced  Low  Emission 
Combustors 


MAX. 

POWER 


Relationship  between  advanced  low  emission 
combustors  and  conventional  combustors  vs. 
1979  EPA  regulations  (revised  from  Verkamp 
et  al.,  1974). 
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Walsh,  1975;  Wairpler  et  al.,  1970)  propulsion  systems.  Here  the 
combustor  of  Fig.  1-1  is  modified  to  include  a fuel  preparation  tube 
located  upstream  of  the  primary  combustion  zone.  The  sole  purpose  of 
the  mixing  tube  is  to  coii5)letely  prevaporize  the  liquid  fuel  and  mix 
the  resulting  fuel  vapor  with  air  prior  to  combustion.  This  not  only 
reduces  thermal  NO^  substantially  (Mellor,  1976;  Verkamp  et  al.,  1974) 
but  also  leads  to  reduced  smoke  and  radiation  levels  as  well  as 
improved  turbine  inlet  tenq^eratures  (Verkanp  et  al.,  1974)  because 
regions  of  near  stoichiometric  equivalence  ratio  and  temperature 
normally  present  in  conventional  combustors  are  eliminated. 

Unfortunately,  the  absence  of  high  local  fuel  concentrations 
and  tenperatures  in  the  primary  combustion  zone  severely  narrows  the 
flame  stabilization  limits  (Verkamp  et  al.,  1974;  Roberts,  P.  et  al., 

1976;  Niedzwiecki,  1975;  Roberts,  R.  et  al.,  1975  and  1976;  Roffe  and 
Ferri,  1976;  Anderson,  1975)  causing  blowoff,  ignition  and  altitude 
relight  problems  over  the  range  of  normal  combustor  operating  condi- 
tions. Furthermore,  since  the  fuel  is  introduced  upstream  of  the 
flameholder,  flame  propagation  and  autoignition  have  been  observed 
in  the  fuel  preparation  tube  (Roffe  and  Ferri,  1975;  Roberts,  P.  et  al., 

1976;  Niedzwiecki,  1975;  Bahr,  1975;  Gleason  et  al.,  1976;  Roffe,  1976; 

Anderson  et  al.,  1975;  Blazowski  and  Bresowar,  1974;  Blazowski  and 
Walsh,  1975;  Wanpler  et  al.,  1976;  Marek  and  Papathakos,  1976; 

Anderson,  1975  and  1976;  Craig,  1976;  Clayton,  1976;  Marek  et  al., 

1977) . Problems  of  inconqjlete  prevaporization  and  premixing  also 
limit  the  effectiveness  of  these  combustors.  Several  investigators 
(Mellor,  1976;  Roffe  and  Ferri,  1975  and  1976;  Roberts,  P.  et  al.,  1976; 

Roberts,  R.  et  al.,  1975  and  1976;  Mosier  and  Roberts,  1974;  Azelbom 
et  al.,  1973;  Blazowski  and  Bresowar,  1974;  Altenkirch  and  Mellor,  1974; 

Marek  and  Papathakos,  1976;  Tacina,  1976)  have  found  it  extremely 
difficiilt  to  prevaporize  and  premix  typical  hydrocarbon  fuels  in  a 
period  of  milliseconds. 

The  present  study  examines  problems  of  blowoff,  flashback, 
inconplete  prevaporization  and  combustion  efficiency  associated  with 
prevaporizing/premixing  combustors  using  a sinplified  axisynmetric  i 

flameholder  configuration  (Chapter  III) . This  combustor  represents  a I 

i 

i 

I 


J 
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coitpromise  between  the  disc  configuration  studied  by  Tuttle  et  al.  (1975, 
1976,  1977)  , the  dunp  ramjet  combustor,  and  the  prevaporizing/premixing 
tube  in  a turbojet.  A fuel  preparation  tube  attached  to  the  front  of 
a disc  permits  at  least  partial  vaporization  of  the  liquid  fuel 
(injected  at  the  upstream  end  on  the  centerline  of  the  tube)  prior  to 
combustion.  The  disc,  mounted  at  the  end  of  the  fuel  preparation  tube, 
stabilizes  the  flame  by  means  of  a bluff  body  recirculation  zone  which 
is  essentially  identical  to  the  flame  stabilization  enployed  in  ramjet 
and  turbojet  afterburners.  Data  obtained  include  variations  in  combustor 
pressure,  inlet  tenperature,  velocity,  geometry  and  equivalence  ratio. 

In  addition,  the  degree  of  prevaporization  is  also  examined  by  varying 
injector  size  and  fuel  type. 

Chapter  II  introduces  the  modeling  approach,  which  will  be 
used  for  the  combustion  efficiency  and  flame  stabilization  data,  and 
shows  examples  of  previous  correlations.  The  model  proposed  by 
Mellor  (1976)  separates  the  combustion  process  into  heterogeneous, 
fluid  mechanic  and  chemical  effects  each  of  which  are  characterized 
by  an  appropriate  time  scale  and  easily  conputed  from  combustor 
inlet  conditions.  Although  these  characteristic  times  are  order  of 
magnitude  estimates,  they  can  be  used  to  identify  the  inportant 
physical  processes  occurring  in  a particular  region  of  the  flow  field. 
Tuttle  et  al.  (1977)  developed  linear  correlations  for  NO^  and  CD 
emissions  using  a sinple  disc  configuration.  Recently,  Mellor  (1977a,b) 
and  Hammond  (1977)  have  been  able  to  extend  these  emission  models  to 
practical  gas  turbine  combustors.  In  addition  to  correlation  of 
existing  experimental  data,  the  characteristic  time  approach  has  been 
used  to  optimize  burner  design  (Mellor,  1977a)  and  predict  quite 
accurately  resulting  NO^  and  CD  emission  levels  (Haimond,  1977) . 

However  the  model  does  not  yet  consider  lean  blowoff  which  was  observed 
under  low  power  conditions  (Hammond,  1977)  or  problems  associated  with 
prevapori zing/premixing  combustors . 

Limited  pressure,  temperature  and  velocity  measurements  were 
performed  (Chapter  IV)  in  the  immediate  vicinity  of  the  fuel  preparation 
tube  to  help  characterize  the  flow  field  for  the  combustion  efficiency 
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(Chapter  V)  and  flame  stabilization  (Chapter  VI)  correlations.  In 
this  study,  characteristic  times  for  lean  blowoff  are  quantified  for 
the  AFOSR  cc.»figuration  as  well  as  data  obtained  by  Ballal  and 
Lefebvre  (1977)  using  a baffle  flameholder.  Following  Zukoski  and 
Marble  (1956) , the  characteristic  time  model  conqjares  a fluid  mechanic 
time  which  is  a function  of  geometry  and  velocity,  to  a fuel  ignition 
delay  time.  Both  times  are  evaluated  in  the  shear  layer  region  between 
the  recirculation  zone  established  by  the  disc  and  the  main  stream. 

The  result  is  a linear  correlation  with  a y- intercept  near  zero  as  the 
theory  predicts.  Analysis  of  the  fluid  mechanic  time  also  predicts 
an  optimum  flameholder  size  for  stability  at  a given  power  setting 
which  is  solely  a function  of  housing  diameter.  In  addition,  the  fluid 
mechanic  time  for  lean  blowoff  is  identical  to  the  one  obtained  by  Tuttle 
et  al.  (1977)  for  minimum  CX)  emissions;  however,  chemical  times  are 
different.  Inconplete  vaporization  is  also  included  in  the  model 
through  a heterogeneous  droplet  evaporation  time. 

Flame  flashback  into  the  fuel  preparation  tube  is  also 
considered  in  Chapter  VI.  Plee  and  Mellor  (1977)  have  reviewed  the 
literature  and  found  that  classical  flashback  has  not  been  observed 
in  non-catalytic  combustors  burning  hydrocarbon  fuels . Instead, 
autoignition  of  the  air  stream,  flame  propagation  through  reversed 
flow  fields  and  preignition  of  the  fuel/air  mixture  in  separated  flow 
regions  of  the  mixing  tube  are  responsible  for  the  reported  phenomena. 
Therefore,  in  non-catalytic  combustors,  classical  flashback  can  be 
avoided  by  eliminating  all  possible  flow  disturbances.  Only  catalytic 
conijustors  operating  at  very  low  reference  velocities  appear  capable 
of  ejqjeriencing  classical  flashback. 

Combustion  efficiency  measurements  have  not  been  con^Jletely 
correlated  to  date  using  the  characteristic  time  model  (Chapter  V) . 
Therefore,  detailed  internal  measurements  are  planned  in  the  future 
(Chapter  VII)  to  help  characterize  regions  of  combustion  inefficiency. 
These  measurements  include  local  species  concentration  and  teii5)erature , 
as  well  as  velocity  magnitude  and  direction  via  LDV. 
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II.  CHARACTERISTIC  TIME  MODEL 
FOR  SPRAY  COMBUSTION 

Modeling  of  liquid  fueled  combustors  is  an  extremely  complicated 
problem  involving  two  phase  turbulent  flow  with  chemical  reaction.  Fluid 
mechanics,  homogeneous  kinetics,  fuel  properties  and  fuel  spray  charac- 
teristics are  all  important  parameters  which  control  the  combustion  process. 
Since  simultaneous  mathematical  modeling  of  the  fluid  mechanic,  chemical 
and  heterogeneous  effects  represents  a formidable  task  at  the  present  time, 
a characteristic  time  model  has  been  developed  as  a convenient  tool  for 
correlating  and  predicting  continuous  combustor*  behavior. 

Other  phenomenological  models  for  high  intensity  gas  turbine 
combustors  exist.  However,  the  characteristic  time  model  appears  capable 
of  including  effects  of  changes  in  not  only  inlet  conditions,  but  also 
combustor  geometry  on  gaseous  emissions,  and  thus  due  to  space  limitations 
will  receive  primary  attention  here.  The  interested  reader  is  referred  to 
recent  reviews  by  Odgers  (1975)  and  Hammond  and  Mellor  (1970)  for  descrip- 
tion of  other  such  models . A complete  model  would  also  predict  smoke 
emissions,  local  heat  flux  to  the  wall,  lean  blowoff  limit,  combustor 
exhaust  plane  temperature  profile  (pattern  factor)  and  so  on. 

The  characteristic  time  model  separates  spray  combustion  into 
three  subprocesses  which  include  1)  vaporization  of  the  liquid  fuel 
droplets,  2)  mixing  of  the  fuel  vapor  with  air,  and  3)  combustion  of  the 
resulting  fuel/air  mixture.  Each  of  these  subprocesses  must  occur  before 
combustion  is  complete.  According  to  Mellor  (1976),  these  subprocesses 
can  be  represented  by  characteristic  times  (Table  2-1)  which  control 
emissions,  performance  and  flame  stabilization. 


"E.g.,  as  in  turbojets  and  ramjets. 
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A.  Characteristic  Times 

The  first  time  listed  in  Table  2-1  is  the  fuel  droplet 

lifetime  which  determines  the  heterogeneous  properties  of  the  combustion 
process.  This  time  is  a function  of  nozzle  characteristics,  fuel  properties, 
fuel  flow  rate  and  ambient  combustor  conditions.  For  a properly  designed 
injector  and  combustor,  should  be  negligible  (Mellor,  1977a, b). 

The  second  and  third  characteristic  times  refer  to  fluid 
mechanic  mixing  processes  occurring  in  the  flow  field,  represents 
the  nixing  process  in  the  region  of  the  fuel  injector  which  is  expected 
to  be  important  in  combustors  utilizing  air  blast  and  air  assist  nozzles 
to  decrease  smoke  emissions.  Only  gaseous  emissions  are  of  concern  in 
the  present  work  (but  see  Dodds  et  al.,  1976);  therefore,  will  be 
omitted  from  further  discussion. 


Table  2-1, 


Time 

Fuel  droplet  lifetime 

Eddy  dissipation  time 
for  injected  fluid 


Eddy  dissipation  time 
in  the  shear  layer 

Fuel  ignition  delay  and 
burning  time 

NO  formation  time 


Characteristic  Times  for  Combustion 
and  Pollutant  Formation  in  Two-phase 
Turbulent  Flow  (from  Mellor,  1976) 


Symbol 

^eb 

^fi 


sJl 


he 


no 


Physical  or  Chemical  Processes 

Droplet  evaporation  and/or  combustion 

Small  scale  turbulent  mixing  near 
the  fuel  injector  in  the  recircu- 
lation zone 

Large  scale  turbulent  mixing  between 
fresh  air  and  the  recirculating 
burned  gas -fuel  mixture 

Homogeneous  combustion  of  the  fuel 
to  CO2 

Homogeneous  kinetics  for  NO  formation 


X characterizes  the  turbulent  mixing  process  in  the  shear 
si 

layer  between  the  fresh  incoming  air  and  the  recirculated  burned 
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gas -fuel  mixture.  For  premixed/prevaporized  combustors,  this  mixing 
time  in  the  shear  layer  would  refer  to  the  turbulent  heat  transfer 
between  the  incoming  mixture  and  the  hot  burned  gases  which  hold  the 
flame. 

The  last  two  times,  and  listed  in  Table  2-1  represent 
homogeneous  kinetic  times  for  the  combustion  of  the  fuel  and  NO  formation. 
These  processes  occur  after  the  fuel  and  oxidizer  have  been  mixed  on  a 
molecular  level  and  are  functions  of  the  local  equivalence  ratio, 
ten^jerature  and  combustor  pressure,  x^^^  may  be  used  to  characterize 
(20  and  HC  oxidation  (combustion  efficiency)  or  HC  ignition  delay  time 
(blowoff  and  flashback).  Exanples  of  previous  work  utilizing  the 
characteristic  time  model  follow. 

B.  Applications  of  the  Model--Disc-in-Duct  Burner 

Under  EPA  sponsorship  Tuttle  et  al.  (1976)  have  extensively 
studied  exhaust  pollutant  emissions  for  the  configuration  shown  in 
Fig.  2-1,  an  axisymmetric  disc- in-duct  burner  with  liquid  fuel 
injection  via  a hollow  cone  spray  nozzle  mounted  at  the  center  of  the 
disc.  Inlet  air  teirperature , pressure,  air  and  fuel  flow  rates,  fuel 
and  nozzle  type  and  disc  geometry  (diameter  and  in  some  cases  teeth  on 
the  outer  circumference)  all  were  varied,  (las  sampling  and  tenperature 
measurements  were  also  made  at  various  points  inside  the  combustor  and 
used  to  develop  the  flow  model  shown  in  Fig.  2-2.  Note  that  the  region 
of  HC  and  CO  quenching  is  separate  from  the  region  of  NO  formation. 

The  correlations  are  developed  from  the  characteristic  times 
as  follows,  assuming  for  the  moment  that  the  droplet  evaporation  time 
x^jj  is  negligible  compared  to  the  shear  layer  mixing  time  x^^^:  for 
a near-stoichiometric  eddy  capable  of  forming  NO,  just  downstream  of 
the  recirculation  zone  shown  in  Fig.  2-2 


eddy 


djNO] 


= 2kj[0][N2] 


2-1 


J 
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where  the  Zeldovich  luechanism  is  assumed  to  apply  and  where  refers  to 
0 + N2  ^ NO  + N.  Invoking  0 ? 5^2  equilibrium  with  equilibrium  constant  K 


dfNOl 

dt 


2-2 


eddy 


T by  definition  is  the  characteristic  kinetic  time  for  NO  formation, 
no  ^ 


proportional  to  K/k^^,  so  that 


[02]''[N2] 

eddy  ^no 


2-3 


Taking  oxygen  and  nitrogen  as  constant  during  the  eddy  lifetime,  as  well 

as  its  tenperature  prior  to  quenching,  Eq.  2-3  may  be  integrated  over  this 

lifetime,  x „ , to  yield 

sx.,no 


[NO] 


eddy  ^sil,no^^no 
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The  total  number  of  NO-forming  eddies  is  assumed  proportional  to  the 
total  fuel  flow  rhj. 


[NO]^  1 '''  T „ /x 

Hotal  f s8.,no  no 
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Thus  the  emissions  index  is  formed 


NO^EI  (gNO^/kg  fuel)  •>. 


2-6 


Tuttle  et  al.  (1976,1977)  used  the  following  estimates  of  the 
characteristic  times 


_ ,.-12  ^135  kcal/mole/RT.  , 
Xno  - 10  e $-1 


msec 


2-7 


where  R is  the  universal  gas  constant  and  Is  the  adiabatic  stoichio- 
metric flame  temperature  (a  function  of  inlet  air  tenperature,  pressure 
and  fuel  type) . 
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^s£,no  '^'^bg 


where  d is  disc  diameter  and  is  an  estimated  velocity  downstream  of 
the  recirculation  zone  in  Fig.  2-2.  Their  correlation  is  shown  in  Fig.  2-3 
for  liquid  propane  injection  (negligible  heterogeneous  effects);  r,  the 
correlation  coefficient,  is  a measure  of  the  goodness  of  the  least  squares 
fit  shown  in  the  figure  and  equals  one  in  the  ideal  case;  a is  the 
standard  deviation. 

As  droplet  evaporation  times  were  increased  (Jet  A fuel),  NO^EI 
decreased,  most  likely  since  fuel  drops  were  thrown  into  the  relatively 
cold  free  stream  flame  zone  in  Fig.  2-2.  Tuttle  et  al.  (1977)  thus  chose 
to  correlate  the  results  in  terms  of 


NO  El  T . /(t  + ax  , ) 

X si^,no  no  eb-^ 


with  a determined  empirically.  The  combined  correlation  for  the  two  fuels 
is  shown  in  Fig.  2-4. 

Droplet  evaporation  times,  were  estimated  as  follows: 
the  basic  physical  relation  which  governs  either  droplet  evaporation 
or  combustion  is  the  "d^  law"  of  Godsave  (1953) 


2 2 

dr  = - Qt 


where  d^  = instantaneous  fuel  droplet  diameter 

djj  = initial  droplet  diameter 

= evaporation  or  combustion  coefficient 
in  the  combustion  region 

t = time . 


By  allowing  d^  to  become  zero,  the  characteristic  time  for  droplet 
evaporation  or  confcustion  derived  from  Eq.  2-10  is 

d ^ 

= B7  • 

d^  is  a function  of  the  fuel  injector  type,  fuel  properties  and  fuel 
flow  rate  vdiile  3^  is  a function  of  fuel  properties  and  ambient  combustor 


sl.no 


6 


989 -T) 
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conditions.  For  a convective  atmosphere,  is  also  related  to  the 
Reynolds  number  (based  on  mean  droplet  size  and  droplet  relative 
velocity)  and  the  Schmidt  number. 

Mellor  (1972)  has  summarized  several  enpirical  relationships  for 
the  dependence  of  3^  upon  Reynolds  number  and  Schmidt  number.  Here  3^  was 
calculated  from  the  Frbssling  equation  (see  Ranz  and  Marshall,  1952) 

3^  = 3q  (1.0  + 0.276  Re^/^Sc^/^)  2-12 

vdiich  corrects  for  forced  convection.  3^  is  defined  as 


3 


o 


8k 

Pjl^p 


Sji  (1  + B) 


with 


2-13 


B = 


2-14 


where  3 = evaporation  coefficient  in  a stagnant 
atmosphere 

k = gaseous  thermal  conductivity 

= density  of  liquid  fuel  at 

Cp  = gaseous  specific  heat  at  constant  pressure 

B = transfer  number  for  droplet  evaporation 

T = tenperature  of  ambient  gas 

= boiling  point  of  liquid  fuel  at  the  local 
pressure 

L = sensible  enthalpy  of  the  liquid  fuel  from 
the  initial  tenperature  to  T.  plus  the 
latent  heat  of  evaporation.  ^ 


Since  practical  liquid  fuel  nozzles  do  not  produce  a uniform 
spray  of  constant  d^,  but  rather  a droplet  distribution,  a special 
mean,  the  Sauter  mean  diameter  (SMD) , is  used  for  the  calculation  of  d 
in  1x4.  2-11.  The  SMD  is  defined  as  the  diameter  of  a droplet  vdiose 
volume  to  surface  ratio  is  the  same  as  the  volume  to  surface  ratio  of 
the  entire  spray.  Thus,  the  SMD  is  a measure  of  the  degree  of  atomization 
and  can  be  conputed  via  empirical  relations  which  are  sumnarized  in  Tuttle 
et  al.  (1975).  For  the  sinplex  pressure  atomizing  nozzles  in  that  study. 
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the  empirical  equation  of  Hunter  et  al.  (1974)  shown  below  was  used  to 
calculate  the  initial  droplet  diameter 

3^  , .O.SV^  p-fOS  ,.20S 


where  SMD  = Sauter  mean  diameter,  micron 
V = viscosity,  centistoke 

r 

p = density,  gm/cm 
Q = volumetric  flow  rate,  lit/hr 
AP  = nozzle  differential  pressure,  atm. 

CD  emissions,  which  contribute  to  combustion  inefficiency,  were 
correlated  in  a similar  fashion.  For  lean  eddies  in  which  CD  (and  HC) 
oxidation  is  quenched  (near  the  disc  edge) 

= -k2[C0][0H]  2-16 

where  k^  refers  to  CO  + OH  ->  CO2  + H.  Again,  the  characteristic  time  for 
CO  oxidation  is  taken  as 


T ^ 1/k- 

CO  2 


but  was  obtained  einDirically  from  the  data  of  Tuttle  et  al.  (1977)  in 
terms  of  inlet  air  temperature 

2200  cal/mole/RT. 

T = e in  msec  . i 

CO 

Arbitrarily  assuming  [OH]  constant  during  the  oxidation  process  yields 
from  Eq.  2-16 


d[CO]  „ dt 

1(351  ^ ^ 

Integrating  over  the  eddy  lifetime  at  the  outer  edge  of  the  shear  layer 
^^sJl,co^ 
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in  O-  -Teo  /t  « 2-20 

[CD]^  si, CO  CO 

[CO]/[CD]^  = e^^^ijCo/Tco  2-21 

vrtiere  [CO]^  is  the  initial  CO  concentration  within  the  eddy,  taken  propor- 
tional to  fuel  flow  rate  so  that 


[CO]/[CO]^  'V  COE  I (g  co/kg  fuel) 


2-22 


Tuttle  et  al.  (1976)  expanded  the  exponential  into  an  equivalent  power 
series  in  and  found  that  a two  term  approximation  was  adequate 


COEI  ■>.  b . c 


2-23 


where  for  the  parameter  ranges  of  interest 


b « c T /t  0 

CO  si, CO 


and  hence 


COEI  (g  co/kg  fuel)  T /t  „ . 2-24 

The  mixing  time  in  the  CO  quenching  region  of  the  flow  (Fig.  2-2) 
is  taken  as  disc  diameter*  divided  by  cold  air  velocity  in  the  annuliis 
surrounding  the  disc, 


T 


si, CO 


2-25 


Corresponding  results  for  C^Hg  fuel  are  shown  in  Fig.  2-5. 

When  Jet  A was  burned  in  the  disc- in-duct  configuration, 
combustion  inefficiency  increased  due  to  the  effective  increase  in  droplet 
lifetime  with  this  heavier  fuel  and  the  consequent  penetration  of  droplets 


*For  discs  with  teeth  a more  conqilicated  scale  factor  was  required;  see 
Tuttle  et  al.  (1976). 
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into  the  cold  free  stream  flame  zone  of  Fig.  2-2.  To  account  for  the 
increase  in  con±)ustion  inefficiency  from  heterogeneous  effects,  the 
correlation  including  the  droplet  lifetime,  took  the  form 

COEI  ■>.  ♦ e 2-26 

where  enpirically  e = 0.24.  COEI  data  for  both  C^Hg  and  Jet  A are 
correlated  with  this  relationship  as  shown  in  Fig.  2-6. 

It  is  inportant  to  realize  these  characteristic  times  are  at 
most  order  of  magnitude  estimates  for  the  various  processes  relevant  to 
the  spray  combustion  flame.  At  best  they  demonstrate  the  proper 
functional  dependence  upon  combustor  design,  cycle  operating  points  and 
fuel  properties.  Thus  no  physical  significance  is  attributed  to  quanti- 
ties such  as  e in  Eq.  2-26,  but  with  enq)irical  determination  of  these 
values,  the  designer  has  available  quantitative  tools  to  coiT5)ute,  in  this 
case,  changes  in  pollutant  emissions  due  to  changes  in  burner  geometry, 
combustor  loading,  and  fuel  type  (see  Tuttle  et  al.,  1976,  1977).  The 
next  section  demonstrates  that  similar  scaling  procedures  are  applicable 
to  conventional  turbojet  liners. 

C.  Applications  of  the  Model --Turbojet  Combustors 

In  the  preceding  sections  the  time  scales  characterizing  the 
dominant  subprocesses  occuring  within  gas  turbine  combustors  were  identi- 
fied from  first  principles.  Emissions  correlations  incorporating  these 
characteristic  times  were  then  derived  and  quantified  for  a confined, 
axisynmetric  disc-stabilized  turbulent  diffusion  flame  burner.  Ifowever, 
the  ultimate  utility  of  the  characteristic  time  model  and  associated 
correlations  lies  in  their  ability  to  correlate  and  hence  predict  exhaust 
plane  emissions  as  a function  of  confcustor  loading,  fuel  type,  and 
burner  geometry,  for  actual  gas  turbine  hardware  exemplified  in  turbojet 
combustors . 

The  disc- in -duct  laboratory  burner  simulates  the  primary  zone  of 
conventional  liners  via  the  existence  of  a central  recirculation  zone,  fuel 
spray-shear  layer  interaction,  and  essentally  identical  but  separate  regions 


23 


of  high  concentrations  of  NO  and  combustion  inefficiency  (see  Fig.  2-2). 

That  is,  in  a conventional  combustor  (see  for  exanple  Mellor,  1976)  NO 
generally  forms  just  downstream  of  the  recirculation  zone(s),  whereas 
HC  and  CO  quench  as  a result  of  air  addition  through  the  liner  walls 
(air  addition  for  the  disc- in-duct  burner  is  limited  to  the  annulus 
surrounding  the  disc).  Therefore,  the  relationships  that  successfully 
correlated  disc-in-duct  emissions  should  extend  to  practical  combustors  . 

In  the  following  the  emissions  correlating  parameters 
developed  above  are  extended  to  two  conventional  jet  stabilized  can  com- 
bustors (Mellor,  1977a, b):  the  General  Motors  GT-309  regenerative 
automotive  burner  (Turunen  and  Collman,  1965;  Wade  and  Cornelius,  1972) 
and  the  Detroit  Diesel  Allison  T-63  helicopter  burner  (Troth  et  al.,  1973,  | 

1974;  Moses,  1975);  and  to  a conventional  swirl  stabilized  annular 

combustor  (Mellor  and  Washam,  1978) , the  Pratt  and  Whitney  JT9D  aircraft  j 

burner  (Nelson,  1976) . | 


! 
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1.  Emissions  Correlating  Parameters  for  Conventional  Combustors 

The  two  conventional  can  combustors  to  v4iich  the  emissions 
correlations  are  extended  are  typified  schematically  in  the  upper 
portion  of  Fig.  2-7  (here  shown  for  the  standard  GT-309  burner).  Of 
significance  in  this  figure  is  the  designation  of  axial  distance  from 
the  injector  face  to  the  respective  discrete  points  of  air  addition 
(i.e.,  slot^  these  lengths  that  will  enable  correlation 

of  geometry  effects  upon  exhaust  emissions.  Tne  annular  JT9D  liner  is 
considerably  more  conplex  but  an  axial  cross  section  taken  through  one 
of  the  injectors  reveals  the  essential  features  depicted  in  the  lower 
portion  of  Fig.  2-7. 

In  dealing  with  actual  gas  turbine  engines  and  conventional 
fuels,  it  is  a reasonable  first  assunption  that  the  injector  provides 
good  atomization  and  that  heterogeneous  effects  are  negligible  (i.e., 

Tgjj  is  effectively  small  conpared  with  x^^^) . Turbulent  mixing  times 
in  the  shear  layer  should  be  coin)uted  as  a ratio  of  local  turbulence 
integral  length  scale  to  turbulent  intensity.  Since  these  quantities 
are  not  available  for  the  liners  of  interest,  mixing  times  are 
evaluated  as  proportional  to  the  ratio  of  a significant  macrodimension 
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to  some  time-mean  velocity  as  for  the  disc- in-duct.  NO^  and  CD  emissions 
are  also  expected  to  originate  from  different  regions  in  these  combustors 
and  hence  t . and  t are  not  necessarily  identical, 

S Xf  j ilO  S X*  j CO 

The  chemical  kinetic  times  are  defined  so  as  to  properly  include 
tenperature  effects  ipon  emissions.  This  tenperature  dependence  is 
realized  in  conjunction  with  an  activation  energy  term,  and  pre- 
exponential factors  are  chosen  to  produce  ratios  of  these  times  with 
respective  mixing  times  of  order  unity. 

Hence,  with  the  above  assunption  of  mixing  controlled  combustion 
and  the  determination  of  the  turbulent  mixing  and  chemical  kinetic  times, 
correlations  of  the  form  found  in  Eq.  2-6  and  2-24  are  applied  to  the 
GT-309,  T-63  and  JT9D  engine  combustors. 

2.  N0„  Correlation  for  Conventional  Combustors 

Lipfert  (1972)  and  Sawyer  et  al.  (1973)  demonstrated  the 
significance  of  the  adiabatic  stoichiometric  flame  temperature  T^_j  on 
NO^  formation  in  several  gas  turbine  engines.  Accordingly,  the  chemical 
kinetic  time  for  NO  formation  x for  the  three  engines  of  interest  is 

based  on  T^_t  and  similar  to  Eq.  2-7. 

(p  X 

The  effect  of  combustor  geometry  upon  NO^  emissions  is 
displayed  in  Fig.  2-8  for  five  different  GT-309  automotive  burners  (Wade 
and  Cornelius,  1972).  The  ability  of  x to  include  inlet  air  temper- 
ature effects  for  a given  liner  is  demonstrated,  but  the  influence  of 
geometry  on  NO^  emissions  remains  (Mellor,  1977a). 

For  the  correlation  to  include  geometry,  determination 

of  the  turbulent  mixing  time  in  the  shear  layer,  x^^^  is  necessary. 
Characteristic  of  the  GT-309  burner  (Fig.  2-7)  are  the  eight  secondary 
air  penetration  holes  which  force  the  shear  layer  to  the  centerline 
and  therefore  promote  closure  of  the  recirculation  zone,  constituting 
a jet  stabilized  flame;  the  air  penetration  tlirough  these  secondary  jets 
quenches  NO  formation  at  this  axial  station.  To  include  the  effect 

X 

of  burner  diameter  it  was  necessary  to  assume  for  the  NO^  macroscale 
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GASIFIER  TURBINE  SPEED  > PERCENT  OF  DESIGN 

Figure  2-8.  NO  (ppm)  versus  percent  of  design  power  for  five 
GT*309  burners  (fran  Wade  and  Cornelius,  1972). 
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v^iere  is  the  axial  distance  from  the  point  of  fuel  injection  to 
the  centerline  of  the  secondary  hole  row  (see  Fig.  2-7)  and  is 

combustor  diameter  at  the  axial  plane  of  the  secondary  jets.  Eq.  2-27 
also  defines  the  JT9D  NO  macroscale  (Mellor  and  Washam,  1978)  and  a 
similar  expression  specifies  that  for  the  T-63  combustor  (Mellor,  1977b). 

The  velocity  in  the  prevalent  NO^  forming  region  of  the 
flowfield  is  taken  as 


m 


a,p  f 
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where  m is  the  primary  zone  air  mass  flow  rate  into  the  contustor 
a ,p 

of  cross  sectional  area  is  combustor  inlet  pressure,  and  M is 

taken  as  the  molecular  weight  of  air.  The  mixing  time  for  NO^  in 
conventional  combustors  is  then 


^sil,no 


= Z 


no 
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NO  El  as  a function  of  x ^ /x  is  shown  in  Fig.  2-9  for  ten 
X sJl,no'  no  ^ 

different  GT-309  burners  and  three  different  fuel  injectors  and  the 
correlation  coefficient  r of  0.932  indicates  fairly  good  correlation. 
Data  collapse  independent  of  geometry  indicates  that  the  correlation 
includes  this  effect  well. 

The  capability  of  x to  include  differing  inlet  conditions 

(represented  by  different  laboratories)  and  also  fuel  type  (with  the 

exception  of  #5  oil*)  upon  NO^EI  for  a fixed  geometry  is  shown  in 

Fig.  2-10  for  the  standard  T-63  liner.  Unfortunately,  the  regenerative 

cycle  was  not  as  effectively  correlated.  Furthermore,  in  that  nozzle 

variations  also  collapsed  to  the  curve  (within  the  scatter)  for  the 

T-63  and  Gr-305  (Fig.  2-9)  coiribustors  confirms  the  initial  assunption 

that  for  these  injectors  distributing  the  listed  fuels  (DF2  was  the  only 

fuel  burned  in  the  GT-309),  the  flame  remains  mixing  controlled 

(i.e.,  X „ » T 1.).  Note  that  the  #5  oil  is  the  only  fuel  of  those 

sx,no  eb 

listed  with  substantial  fuel-bound  nitrogen.  Geometry  effects  for 


*It  is  the  convention  here  to  flag  data  that  are  excluded  from  the 
least  squares  fit. 
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conventional  T-63  combustors  were  correlated  with  a result  similar  to 

that  shown  in  Fig.  2-9  (see  Mellor,  1977b).  NO  El  data  from  the  JT9D 

Mod  II  combustor  geometry  also  correlated  quite  well  with  t /t  . 

sj^,no  no 

Data  supplied  indicated  that  t also  includes  the  effect  of  ambient 

no 

humidity  iqx)n  NO^  emissions  (Mellor  and  Washam,  1978) . 

3.  ro  Correlation  for  Conventional  Combustors 

In  a gas  turbine  combustor  it  is  the  discrete  addition  of 
cooler  inlet  air  through  penetration  holes  and  film  cooling  slots 
that  quenches  the  CO  oxidation  reaction  in  the  outer  regions  of  the 
spray  flame.  The  aero thermodynamic  parameters  in  this  region  of  the 
liner  flowfield  then  are  incorporated  in  the  CO  emissions  correlations 
for  engines. 

The  kinetic  time  for  CO  oxidation  t in  liners  is  a function 

CO 

of  the  adiabatic  burned  gas  tenperature  at  the  overall  equivalence 
ratio,  Tj^.  Results  reported  for  the  GT-309  combustor  series  gave 
(Mellor,  1977a) 


= 10  ^ exp  (7700  cal/mole/RTj^)  msec  2-30 

whereas  those  for  the  T-63  configurations  required  the  mean  of  inlet 
air  and  burned  gas  tanperatures  T in  the  CO  kinetic  time  (Mellor,  1977b) 
such  that 


= 10  ^ exp (10, 760  cal/mole/RT)  msec  . 2-31 

Since  the  burned  gas  tenperature  is  primarily  a function  of 
inlet  air  ten?)erature , equivalence  ratio  and  fuel  type,  these  effects 
upon  CO  emissions  are  included  through  Evidence  of  this  is  provided 

in  Fig.  2-11  where  COEI  data  corresponding  to  different  operating  condi- 
tions and  four  different  fuels  for  the  standard  T-63  liner  are  plotted 
versus  x /x„.  . With  the  exception  of  #2  Diesel  fuel  and  #5  oil 

(flagged  in  this  figure) , data  for  JP4  and  Jet  A collapsed  to  the 
mixing  controlled  correlation  independent  of  injection  scheme  signifying 
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that  heterogeneous  effects  for  these  cases  are  negligible.  Different 
inlet  conditions  (the  two  laboratories)  as  well  as  both  non-  and 
regenerative  cycles  are  also  collapsed  by  the  correlation.  The  increase 
in  combustion  inefficiency  with  #2  Diesel  fuel  and  #5  oil  is  attributable 
to  the  significant  droplet  effect  associated  with  these  relatively  heavier 
fuels  (Mellor,  1977b). 

Determination  of  will  permit  inclusion  of  geometry 

effects  ipon  CO  emissions.  The  CO  macroscale  for  liners  £ is  identical 

CO 

in  form  to  that  for  NO^EI  given  previously,  where  axial  and  radial 
influences  are  included  through 


^co  ^quench  ^ '^comb  ‘ 

However,  the  CO  macroscale  exhibits  somewhat  more  complicated  behavior 

than  that  for  NO^.  Specifically,  the  GT-309  combustor  series  experienced 

an  axial  quench  shift  for  CO  with  power  setting  (i.e.,  decreased 

with  decreasing  engine  power).  For  high  power  (70  to  100%  of  design) 

CO  quenched  in  the  vicinity  of  the  secondary  holes.  Hence,  for  high 

power,  £^^  was  based  upon  At  low  power  (50  to  60%  of  design)  CO 

quenched  in  either  the  region  of  the  dome  holes,  film  cooling  slot  or 

the  primary  holes  depending  upon  combustor  configuration  (geometry) . 

£„  for  low  power  conditions,  then,  was  based  on  £,  , £j..,  , or  £ . 

CO  ’ ’ dome  film  pri 

(Mellor,  1977a).  Several  T-63  configurations  also  demonstrated  this 
axial  CO  quench  shift  with  power. 

The  turbulent  intensity  relevant  to  CO  quenching  in  the  GT-309 
burners  was  taken  as  proportional  to  the  mean  air  injection  velocity 
through  the  holes  and  slots  in  the  liner  (Mellor,  1977a) 

''a.holes  = ” \oles  ^-3 

where  Aj^g^gg  is  tl;e  total  open  area  of  the  particular  combustor.  In  the 
T-63  series  combustors  turbulent  intensity  scaled  with  combustor 
reference  velocity  instead  of  air  injection  velocity 
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That  is,  for  the  T-63  burners,  no  dependence  of  COEI  on  liner  open  area 
was  found  (Mellor,  1977b).  Geometry  effects  upon  CX)  emissions  from 
these  two  liners  are  thus  included  by 
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T-63  COEI data  are  displayed  in  Fig.  2-12  where  (excluding  the  flagged 
advanced  burner  data  and  by  accounting  for  the  CO  quench  shift)  all 
geometry  variations  listed  are  observed  to  collapse  rather  well  to  the 
curve.  GT-309  COEI  data  also  experienced  excellent  geometry  correla- 
tion by  allowance  for  the  quench  shift. 


4.  Predictions  and  Validation  of  tiie  Model 

With  relatively  high  statistical  levels,  NO^  and  CO  emissions 
from  mixing  controlled  flames  within  the  GT-309  combustor  series  have 
been  correlated  and  hence  emissions  predictions  for  this  family  of  geo- 
metrically similar  burners  are  possible.  Specifically,  since  these 
correlations  include  variations  in  combustor  loading  as  well  as  liner 
geometry,  curves  particularly  relevant  to  the  combustor  designer  are 
produced  (each  at  a given  power  setting)  displaying  the  effect  of 
variations  in  liner  configuration  upon  exhaust  emissions. 

Fig.  2-13  is  a plot  of  both  NO^EI  and  COEI  versus  the 
non- dimens ionali zed  geometry  parameter  90%  power  point. 

The  discrete  symbols  are  actual  data  taken  on  the  listed  burners  at 
General  Motors.  In  testing  configurational  variations  on  the  baseline 
and  decreased  diameter  burners,  two  diameters  were  effectively  examined. 
This  graph  demonstrates  the  significant  effect  of  combustor  diameter 
upon  high  power  NO^:  burner  diameter  can  be  reduced  to  achieve  a sub- 
stantial decrease  in  high  power  NO^  without  appreciable  sacrifice  in 
high  power  CO.  Furthermore  ^ ®3n  be  diminished  to  approximately 

0.3  to  effect  high  power  NO^  reduction  before  incurring  a substantial 
increase  in  high  power  CO. 
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Figure  2-12.  COFI  versus  Tco/t^o  for  several  T-63  burners 
(Mel lor,  1977b).  ’ 
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Similar  graphs  for  given  power  points  are  available  from 
the  GT-309  correlations.  A plot  for  60%  power  is  also  examined  in 
Mellor  (1977a)  and  demonstrates  the  requirement  of  large  combustor 
diameters  and  delayed  air  addition  for  minimizing  low  power  CO;  both 
trends  are  opposite  to  those  necessary  for  reducing  high  power  NO^. 

Hence  to  minimize  both  NO^  and  CO  emissions  from  a given  fixed  geometry 
liner  will  result  in  a geometry  (and  therefore  emissions)  con5)romise 
at  best.  Moreover,  these  results  strongly  suggest  that  the  optimal  low 
emissions  combustor  will  require  variable  geometry. 

Results  from  these  plots  were  then  incorporated  to  design  a 
sanple  optimum  (i.e.,  low  NO^)  Gr-309  burner.  This  sanple  design  was 
based  on  the  standard  GT-309  combustor  so  that  a con^arison  existed 
between  this  optimized  design,  the  known  NO^  benefits  from  the  decreased 
diameter  version,  and  the  unaltered  standard  burner. 

The  essential  feature  of  the  sanple  optimum  is  the  deletion  of 
the  primary  hole  row  to  accommodate  upstream  placement  of  the  secondary 
jets  (to  minimize  high  power  NO  ).  These  relocated  secondary  holes 
are  enlarged  slightly  so  that  standard  burner  mean  air  injection 
velocity  is  maintained. 

The  results  predicted  that  the  sample  optimum  design  should 
exhibit  only  slightly  higher  NO^  levels  over  the  decreased  diameter 
but  without  the  high  levels  of  CO  produced  by  the  latter.  These 
results  also  suggest  that  at  best  a 50%  decrease  in  high  power  NO 
is  expected  from  conventional  (non-prevaporizing/premixing)  gas  turbine 
liners  (for  details  of  sanqjle  optimum  design  and  results  consult  Mellor, 
1977a) . 

To  verify  the  emissions  predictions  from  the  sample  optimum 
burner  three  GT-309  configurations  were  tested  at  General  Motors 
Research  Laboratories  (Hannnond,  1977) . These  were  the  standard  burner  and 
and  two  modifications:  Mod  A,  the  optimum  low  NO^  burner  designed  above 
and  Mod  B,  essentially  Mod  A but  with  four  additional  secondary  jets  of 
the  same  diameter  as  those  in  Mod  A provided  to  examine  the  predicted 
effect  of  mean  air  injection  velocity  i^pon  OOEI. 
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Rig  test  results  for  these  burners  are  shown  first  for  NO  El  in  i 

X j 

Fig.  2-14:  the  predictions  for  the  low  NO^  burners  were  accurate  to 

within  one  standard  deviation  of  test  results.  Hammond  (1977)  attributes 

the  underpredictions  of  NO^  at  high  power  for  the  standard  burner  to  j 

incomplete  NO  quenching  by  the  secondary  jets,  inp lying  that  the  axial  i 

^ 1 

length  scale  for  these  points  is  effectively  greater  than  j 

COEI  emissions  from  these  burners  are  displayed  in  Fig.  2-15.  1 

Agreement  between  predictions  and  experimental  data  is  not  as  good  as  j 

1 

with  NO^.  In  general , predictions  do  not  exist  within  one  standard 
deviation  of  test  data.  Above  70%  gasifier  speed  predicted  and  test  j 

points  exhibit  disparate  trends.  Upon  examination  of  the  actual  test  1 

trends  above  70%  power  (particularly  emphasized  in  combustors  A and  B) , i 

the  data  suggest  that  the  axial  CO  length  scale  increased  continuously  j 

with  power  (i.e.,  it  does  not  attain  an  upper  limit  of  at  and  j 

above  70%  power).  In  fact,  according  to  Hammond  (1977),  CO  quenching  ' 

occurs  at  all  discrete  locations,  the  position  being  dictated  by  | 

power  setting.  j 

Also  indicated  in  Fig.  2-15  is  the  unstable  operation  j 

experienced  by  combustors  A and  B (as  would  be  expected  intuitively  | 

from  the  large  and  early  air  addition  supplied  through  the  secondary 
jets).  Mod  B was  marginally  stable  at  60%  gasifier  speed  and  blowout  | 

occurred  for  both  burners  at  50%  speed.  The  characteristic  time  model  ] 

correlation  did  not  predict  the  occurence  of  blowoff  and  hence  refining  j 
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Figure  2-14.  Conparison  of  predicted  and  experimental  El (NO  ) versus 
gasifier  speed  (from  Hanriond,  1977) . ^ 
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Figure  2-15.  Comparison  of  predicted  and  experimental  El (GO)  versus  gasifier 
speed  (from  Hammond,  1977) . 
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III.  DESIGN  AND  OPERATION  OF  THE  EXPERIMENTAL 
COMBUSTOR  FACILITY 

The  present  combustion  facility  is  equipped  for  testing 
conventional  gas  turbine  combustors  and  simplified  laboratory  burners  such 
as  the  tube-and-disc  configuration  of  this  study.  Instrumentation  allows 
continuous  monitoring  of  all  pressures,  temperatures  and  flow  rates  as 
well  as  the  important  product  species  concentrations. 

Sections  A and  B of  this  Chapter  review  the  basic  system 
hardware  which  has  been  used  extensively  for  other  combustor  programs. 
Details  of  the  facility  are  found  in  descriptions  by  Mellor  et  al.  (1972), 
Tuttle  et  al.  (1973),  Shisler  et  al.  (1974),  Tuttle  et  al.  (1975),  Dodds 
et  al.  (1975)  and  Tuttle  et  al.  (1976).  The  remaining  four  sections 
concentrate  on  the  design  of  the  Air  Force  experimental  apparatus  and  the 
selection  of  run  conditions. 

A.  Facility  Design 

Air  flow  for  the  burner  is  supplied  by  a compressed  air 
blowdown  system  with  a storage  volume  of  approximately  85  m^.  A system 
pressure  of  150  atm  permits  air  flow  rates  of  3 kg/sec  for  40  min;  smaller 
flow  rates  are  generally  used,  however,  which  permit  operation  for  corre- 
spondingly longer  periods  of  time.  The  system  is  capable  of  flow  rates  of 
4.5  kg/sec.  The  air  flow  rate  is  monitored  via  a 5.40  cm  orifice  type 
flowmeter  and  regulated  by  a pneumatically  controlled  air  throttle  valve. 

A schematic  of  the  experimental  apparatus,  comprised  of  two 
combustors  in  series,  is  shown  in  Fig.  3-1.  Commercial  grade  liquid 
propane  is  burned  in  the  slave  combustor  (a  standard  Allison  T-56  liner 
with  a dual  orifice  pressure  atomizing  nozzle)  to  provide  heated  vitiated 
air  (of  temperatures  typical  of  current  air  breathing  engine  operation) 
to  the  test  combustor.  Turbine  flowmeters  measure  the  fuel  flow  to  both 
the  preheat  and  test  combustors.  A remotely  operated  back  pressure 
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Schanatic  of  the  experimental  apparatus. 
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valve  consisting  of  a movable,  conical  centerbody  controls  the  pressure 
in  the  test  section.  Before  release  to  the  atmosphere,  the  exhaust  flow 
noise  level  is  reduced  by  a coirmercial  gas  turbine  silencer. 

The  test  combustor  inlet  ten^erature  is  monitored  with  a 
thermocoi^le  rake  mounted  across  the  diffuser  section.  This  rake  consists 
of  eight  unshielded  chromel-alumel  thermocouples  positioned  in  a stain- 
less steel  support  such  that  they  provide  an  area-averaged  temperature. 

The  thermocouple  junctions  are  insulated  from  the  support  by  ceramic 
connectors  while  the  thermocouple  wires  inside  the  support  are  protected 
by  an  AI2O2  ceramic  filling.  The  fuel  flow  rate  to  the  slave  burner 
controls  the  test  combustor  inlet  temperature. 

Immediately  downstream  of  the  temperature  rake,  gas  samples 
are  extracted  from  the  preheat  combustor  via  an  uncooled  gas  sampling 
rake  (cross -shaped)  constructed  of  stainless  steel  tubing.  The  rake 
has  a total  of  sixteen  collection  holes  spaced  so  as  to  yield  an 
area- averaged  sample. 

A water  cooled  stainless  steel  rake  (described  by  Shisler  et 
al.,  1974)  provides  area-averaged  gas  samples  from  the  exit  of  the  test 
combustor  to  the  gas  analysis  system.  This  rake  is  maintained  at  32.2  cm 
from  the  injector  plane  and  is  inserted  through  the  elbow  section. 

Sampling  is  accomplished  through  sixteen  sampling  holes,  four  located 
along  each  of  the  four  spokes  which  extend  from  a central  hub. 

A water  cooled  single  point  stainless  steel  gas  sampling  probe 
(also  described  by  Shisler  et  al . , 1974)  may  also  be  inserted  in  the 
elbow  section.  The  single  point  probe  is  constructed  of  three  concentric 
stainless  steel  tubes.  Cooling  water  flows  in  and  out  of  the  two  outer- 
most passages  vdiile  the  gas  sample  flows  through  the  inner  passage  to 
the  gas  analysis  system;  the  collected  sanple  is  quenched  in  the  probe 
by  a converging  nozzle  and  the  cooling  water.  The  probe  tip  is  offset 
from  the  main  body  so  that  it  sweeps  out  an  arc  as  it  is  rotated.  The 
probe  can  be  displaced  axially  within  the  test  combustor  allowing  gas 
sample  extraction  from  any  point  along  the  arc  at  any  axial  plane.  A 
platinum/platinum-101  rhodium  thermocouple  mounted  on  the  side  of  the 
probe  provides  an  estimate  of  local  gas  temperatures. 
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This  system  allows  direct  simulation  of  run  conditions  typical 
of  realistic  air  breathing  engines  with  the  exception  of  the  vitiated 
inlet  air.  Table  3-1  summarizes  the  capability  of  this  combustor  facility. 
Note  that  over  the  range  of  operating  conditions  described  in  Table  3-1, 
it  is  possible  to  vary  pressure,  air  flow  rate,  equivalence  ratio,  inlet 
tein)erature  and  fuel  type  independently. 


Table  3-1.  Test  Combustor  Operating  Ranges 


Inlet  temperature 
Combustor  pressure 
Air  flow  rate 
Fuel  flow  rate 
Fuel 


273-950K 

slightly  > 1 atm- 10  atm 
0.5-4. 5 kg/sec 
0.15-2.1  Jl/min 

propane.  Jet  A or  any  typical 
jet  or  Diesel  fuel  (fuel  preheat 
tenperatures  up  to  370K) 


B.  Gas  Analysis  System 

Gas  sanples  are  drawn  from  the  preheat  and/or  test  combustor 
for  analysis  (see  Fig.  3-1).  The  sanple  flow  proceeds  to  the  instru- 
mentation first  through  a fiberglass  filter  for  carbon  particulate 
removal  (98.7%  collection  efficiency  for  0.05  micron  particles)  and  then 
through  8 m of  stainless  steel  sample  line  maintained  at  425K  (Vaught 
et  al.,  1971)  by  thermostatically  controlled  heat  tapes.  Sanples  are 
analyzed  on  a continuous  basis  for  concentrations  of  carbon  monoxide 
(CO),  unbumed  hydrocarbons  (HC) , nitric  oxide  and  nitrogen  dioxide . 

00  and  HC  are  monitored  with  a nondispersive  infrared  analyzer  (Beckman 
model  315A)  and  a flame  ionization  detector  (Beckman  model  402) , respec- 
tively. A chemiluminescent  NO  analyzer  (similar  to  that  of  Fontijn 
et  al.,  1970)  equipped  with  a NO2  stainless  steel  converter  allows 
alternate  monitoring  of  NO  and  total  NO^  (NO  + NO2) . NO2  concentrations 
are  measured  directly  by  a nondispersive  ultraviolet  analyzer  (Beckman 
model  255B) . 
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Each  gas  analysis  instrument  produces  an  output  signal  as  a 
meter  deflection  which  is  observed  visually  and  manually  recorded.  In 
addition,  a continuous  and  permanent  record  of  output  signals  is 
available  from  strip  chart  recorders  connected  to  each  instrument  (similar 
records  are  obtained  for  temperature  monitoring  theimocouples) • 

C.  Combustor  Design 

The  test  combustor  housing  shown  in  Fig.  3-2  is  constructed 
of  stainless  steel  pipe  of  inner  diameter  14.63  cm  and  wall  thickness  of 
1.1  cm.  Five  quartz  windows  (6.35  mm  thick,  3.45  cm  diameter  and 
6.35  cm  apart)  exist  on  both  sides  of  the  housing  to  allow  visual  access 
to  the  test  section.  The  reducing  section  is  wrapped  with  a water  cooling 
jacket.  High  pressure  water  is  injected  downstream  of  all  sanpling 
devices  at  the  location  of  minimum  area  to  protect  the  combustor  facility 
from  high  tenperature  exhaust  gases. 

The  flameholder  geometry  as  positioned  within  the  test  section 
housing  is  shown  schematically  in  Fig.  3-3.  This  configuration  was 
selected  for  the  experiments  and  represents  a conpromise  between  the 
disc-and-duct  configuration  studied  by  Tuttle  et  al.  (1976)  for  EPA 
(see  Chapter  II),  the  dump  ramjet  burner  of  Stull  et  al.  (1974)  and  the 
prevaporizing/premixing  fuel  feed  tube  in  a turbojet.  Liquid  fuel  (Jet  A 
or  propane)  is  injected  into  the  center  of  the  tube  via  a simplex  pressure 
atomizing  nozzle.  The  tube  (length  typical  of  a ramjet  combustor) 
provides  a region  wherein  the  fuel  can  be  partially  prevaporized  and 
premixed  before  combustion.  The  disc  (mounted  in  the  center  of  the 
third  window)  serves  to  stabilize  the  flame  by  means  of  a bluff  body 
recirculation  zone  which  is  essentially  identical  to  the  flame  stabili- 
zation technique  used  in  ramjet  and  turbojet  afterburners. 

To  understand  the  inqxtrtant  effect  of  flameholder  geometry 
upon  both  combustion  efficiency  and  flame  stabilization,  five  geometry 
variations  of  the  basic  tube-and-disc  configuration  are  utilized.  The 
specific  dimensions  of  these  variations  are  delineated  in  Table  3-2 
where  the  variables  correspond  to  those  depicted  in  Fig.  3-3. 
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Table  3-2.  AFOSR  Burner  Geometries 
(See  Fig.  3-3) 


D(cm) 

d(cm) 

L(cm) 

D-d(cm) 

A 

12.7 

4.5 

14.5 

8.2 

B 

10.8 

4.5 

14.5 

6.3 

C 

12.7 

4.5 

8.1 

8.2 

D 

11.7 

3.5 

14.5 

8.2 

E 

12.7 

3.5 

14.5 

9.2 

Geometry  A 's  taken  as  the  reference  or  baseline  burner  and 
was  utilized  for  the  selected  internal  pitot  tube  measurements. 

A streamlined  version  of  the  tube-and-disc  flameholder  is 
shown  in  Fig.  3-4  and  is  essentially  the  baseline  geometry  with  the 
addition  of  a shroud.  This  shrouded  combustor  was  incorporated  primarily 
to  obtain  internal  pressure  and  temperature  measurements  (with  limited 
combustion  efficiency,  blowoff  and  pitot  tube  data)  for  conpafison  with 
the  standard  configuration.  The  shrouded  combustor  was  designed  in 
accordance  with  the  suggestions  of  Edelman  (1976)  to  aid  in  computer 
modeling  of  rairijet  combustors. 

Either  tube-and-disc  assembly  is  mounted  within  the  housing 
by  thin  support  struts  made  of  3 mm  thick  stainless  steel  and  oriented 
so  as  to  minimize  flow  disturbances  in  the  air  flow  stream.  Both  the 
preheat  and  test  combustors  are  ignited  with  CH^/O^  torch  igniters. 

The  fuel  injectors  used  in  this  investigation  are  conventional 
stainless  steel  sinplex  pressure  atomizing  nozzles  with  volumetric  flow 
ratings  of  16  and  30  gal/hr  (water)  at  125  psid  with  a nominal  hollow 
spray  cone  angle  of  90°. 

D.  Combustion  Efficiency  and  Flame  Stabilization  Studies 
Ideally,  operating  conditions  are  chosen  which  would  typify 
advanced  combustors  and  which  would  vary  all  of  the  important  charac- 
teristic times  independently.  However,  this  is  very  difficult  to  do 
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initially  because  the  dependence  of  combustion  efficiency  and  flame 
stabilization  upon  the  characteristic  times  has  not  been  determined  for 
the  tube-and-disc  burner.  Furthermore,  exact  duplication  of  pressure, 
tenperature , reference  velocity  and  overall  equivalence  ratio  for  a 
given  combustor  based  on  another  of  different  geometry  may  result  in 
operating  conditions  outside  of  the  lean  limit. 

As  an  alternative,  the  loading  parameter,  1/0,  defined  by 
Lefebvre  (1966)  is  available  to  simulate  operating  conditions  in  com- 
bustors of  different  geometry.  1/0  vdiich  is  based  on  turbulent  burning 
velocity  theory  is  defined  as 


1/0 


m 

a 

;or 


exp(Tj„/300) 


3-1 


where  P = combustor  pressure,  atm 
T.  = inlet  air  tenperature,  K 

m = air  mass  flow  rate,  kg/sec 

® 2 
\e£  ~ combustor  cross  sectional  area,  m 

= maximum  combustor  width,  m. 

For  this  study,  values  of  l/e  between  0.46  and  3.91  are  used  to  simulate 
the  coitbustor  operating  conditions  of  Blazowski  and  Bresowar  (1974) . 

Having  now  defined  the  combustor  operating  range  in  terms  of 
1/0,  it  is  necessary  to  choose  a base  point  determined  by  the  physical 
limitations  of  the  system  and  the  desire  to  achieve  a prevaporized 
mixture  with  propane  fuel  at  the  exit  plane  of  the  fuel  preparation  tube. 
Consideration  of  the  inportant  characteristic  times  in  the  system  provides 
a means  of  selecting  the  base  point  and  other  operating  parameters. 

Since  conplete  prevaporization  is  not  always  achieved  in 
these  combustors  (Mellor,  1976;  Altenkirch  and  Mellor,  1974;  Roffe  and 
Ferri,  1976),  it  is  desirable  to  estimate  the  droplet  evaporation  time, 
in  the  fuel  preparation  tube.  To  facilitate  pre vaporization,  d^ 

(Eq.  2-11)  is  minimized  by  selecting  a low  fuel  flow,  high  injection 
pressure  nozzle  which  will  provide  enough  flow  to  stabilize  the  flame 
(see  Eq.  2-15).  High  inlet  air  tenperatures  and  low  air  velocities  are 
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also  inportant.  Exact  tube  velocities  to  estimate  droplet  residence 
time  in  the  tube  requires  the  aid  of  experimental  data  (see  pitot  tube 
measurements).  Hence,  combustor  operating  ranges  are  chosen  which 
qualitatively  meet  the  conditions  described  above. 

To  assist  prevaporization,  the  base  point  to  produce  the 
smallest  possible  initial  droplet  diameter  is  selected  as  shown  in 
Table  3-3. 

Table  3-3.  Base  Operating  Point 

Inlet  tenperature  = BOOK  Overall  equivalence  ratio  =0.29 

Combustor  pressure  = 5 atm  Nozzle  = 16-90,  SMD  ==  95  ym  (Jet  A) 

Air  flow  rate  = 1 kg/ sec  Conventional  fuels  = Jet  A 

and  Propane 

From  this  point,  equivalence  ratio,  pressure  and  inlet  temperature  are 
varied  independently  to  determine  the  functional  dependence  of  the  char- 
acteristic times  and  (velocity  variations  were  accomplished  by 
changes  in  pressure  rather  than  air  flow  rate) . is  varied  by 

selecting  a nozzle  (30-90)  which  produces  a much  larger  initial  droplet 
size  ('\'50  ym  larger)  and  hence  changes  the  prevaporization  conditions  in 
the  fuel  preparation  tube. 

Exit  plane  emission  data  relevant  to  combustion  efficiency 
were  collected  with  the  water  cooled,  area-averaging  rake.  Combustion 
efficiency  is  then  quantified  (see  Eq.  5-1)  for  each  run  point. 

Flame  stabilization  data  are  obtained  in  a somewhat  different 
manner.  First,  the  test  cont)ustor  pressure,  inlet  temperature,  and 
air  and  fuel  flow  rates  are  adjusted  to  establish  a well  stabilized 
flame.  After  this  condition  is  maintained  for  a short  time,  the  fuel 
flow  rate  is  slowly  reduced  to  approach  and  ultimately  reach  the  lean 
blowoff  limit  or  increased  in  search  of  flashback.  Experimentally, 
blowoff  is  the  point  where  the  flame  becomes  visually  extinct  whereas 
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flashback  is  the  point  at  which  a thennocouple  mounted  in  the  center  of 
the  tube  near  the  wall  exceeds  a tenperature  of  1200K.  If  flashback 
occurs,  the  fuel  flow  is  automatically  terminated  and  remains  so  until 
the  fuel  supply  system  is  manually  restarted. 

Table  3-4  summarizes  the  range  of  burner  conditions  used  for  the 
combustion  efficiency  and  flame  stabilization  studies.  Selected 
combustion  efficiency  and  lean  blowoff  data  were  also  taken  with  the 
shrouded  combustor  under  the  conditions  listed  in  Table  3-5.  Comparison 
with  the  unshrouded  tube-and-disc  geometry  is  presented  in  Cliapters  V and  VI, 
The  matrix  shown  in  Fig.  3-5  summarizes  the  burner  operating  conditions 
used  in  this  study  as  functions  of  equivalence  ratio  and  1/0. 

E.  Pressure,  Temperature  and  Velocity  Studies 
In  addition  to  the  combustion  studies  discussed  in  the  previous 
section  pressures,  temperatures  and  velocities  inside  and  at  the  exit  of 
the  fuel  preparation  tube  are  of  interest, 

1.  Static  Pressure  and  Temperature  Measurements 

The  fluid  mechanics  of  the  approach  flow  in  the  tube  is  especially 
important  for  the  flashback  study  since  classical  flashback  is  a boundary 
layer  phenomenon  which  may  occur  near  the  tube  exit.  Also,  hot  and  cold 
flow  pressure  measurements  will  show  any  effects  of  combustion  on  the 
velocity  at  the  disc.  To  examine  these  effects,  a differential  pressure 
is  measured  between  stations  0 and  8 (Fig.  3-6)  in  both  hot  and  cold 
flow.  If  the  pressure  drop  remains  constant  for  both  flows  it  is  assumed 
that  the  entrance  flow  prior  to  combustion  is  identical  to  that  in  cold 
flow. 

For  the  remaining  measurements,  the  shrouded  tube-and-disc 
assembly  (Fig.  3-4)  was  instrumented  as  shown  in  Fig.  3-6.  Of  the  four 
static  pressure  measurements,  one  (point  1)  is  an  absolute  pressure  while 
the  other  three  are  differential  pressures  with  point  1 as  reference. 
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Table  3-4.  Combustion  Efficiency  and  Flame  Stabilization 
Burner  Operating  Ranges 

1.0  £ Air  Flow  Rate  £ 1.5  kg/sec 
Fuels  = Jet  A,  Liquid  C^Hg 

2.3  £ Pressure  £ 8.0  atm 

600  £ Inlet  Air  Temperature  £ 800K 

Lean  Blowoff  Limit  £ Overall  Equivalence  Ratio  £ 0.39 

Lean  Blowoff  Limit < Tube  Equivalence  Ratio  £ 2.80 

12.0  £ Reference  Velocity  £ 75.0  m/sec 

31.0  £ Air  Velocity  at  Tube  Exit  £ 190.0  m/sec 

80  £ Initial  Fuel  Spray  Sauter  Mean  Diameter*  £ 160  ym 
0.7  < Mean  Air  Residence  Time  in  Tube  < 6.1  msec 


Table  3-5.  Burner  Operating  Conditions  for  Combustion 
Efficiency  and  Flame  Stabilization  Using  a 
Shrouded  Combustor  (Fuel  = Jet  A;  #16-90  Nozzle) 


Lefe 

Loac 

1/0 

‘bvre 

ling 

1 (1966) 
Parameter 
kg 

Air  Flow 
Rate 

m (kg/sec) 

d 

Inlet  Air  Combustor  Condition 

Temperature  Pressure 

T^  (K)  P (atm) 

sec 

^ iT75  2.75 
atm  m 

1.55 

1.0 

800 

4 

Eff  <1)  = 0.29 

1.05 

1.0 

800 

5 

Eff  = 0.29 

.46 

1.0 

800 

8 

Eff  4)  =0.29 

1.55 

1.0 

800 

4 

b.o. 

.46 

1.0 

800 

8 

b.o. 

K — 

Under  all  test  conditions  except  perhaps  at  8 atm,  liquid  C3Hg  is 
expected  to  flash  vaporize  in  the  fuel  preparation  tube. 


Figure  3-5.  Sunmary  of  burner  operating  conditions. 
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Pressure  taps  are  constructed  of  1/8"  stainless  steel  tubing  and 
securely  welded  to  the  tube  and  the  disc.  Model  A- 10  Sensotec 
0-250  psig  and  ±2  psid  pressure  transducers  monitor  the  absolute  and 
differential  pressures  respectively;  tl\e  combustor  differential  pressure 
is  measured  with  a mercury  manometer. 

Static  temperatures  are  measured  with  Omegaclad  thermocouple 
wire  which  is  constructed  of  1/8"  stainless  steel  sheathing  with  magne- 
sium oxide  insulation  and  24  gage  chromel-alumel  wire.  The  thermocouple 
sheaths  are  also  securely  welded  to  the  disc  and  tube  and  are  mounted  such 
that  the  thermocouple  beads  are  flush  with  the  combustor  surface. 

Table  3-6  summarizes  the  operating  conditions  for  the  instrumented 
flameholder.  Tliese  include  measurements  at  two  different  reference  veloc- 
ities under  conditions  of  cold  flow,  hot  flow  (preheat  combustor  only) 
and  hot  reacting  flow  (both  combustors  operating) . 


Table  3-6.  Burner  Operating  Conditions  for  Static 
Pressure  and  Ten^Derature  Measurements 
with  the  Shrouded  Combustor 


m (kg/sec) 

a 

P (atm) 

T.  (K) 

Vref  Cm/sec) 

Combustion 

1.0 

1.5 

300 

33.8 

none 

1.0 

4.0 

800 

33.8 

preheat  only 

1.0 

4.0 

800 

33.8 

both  combustors 

1.0 

3.0 

300 

16.9 

none 

1.0 

8.0 

800 

16.9 

preheat  only 

1.0 

8.0 

800 

16.9 

both  combustors 

2.  Pitot  Tube  Measurements 

Calculation  of  the  droplet  evaporation  time  in  the  fuel 
preparation  tube  is  possible  only  with  knowledge  of  the  air  flow  split 
or  velocity  at  the  tube  exit.  For  this  reason,  pitot  tube  measurements 
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were  taken  on  the  baseline  tube -and- disc  combustor  under  cold  flow 
conditions  at  the  exit  of  the  tube  and  in  the  annular  area  between  the 
combustor  housing  and  the  disc. 

The  stagnation  pressure  is  measured  with  the  pitot  tube  as 
illustrated  in  Fig.  3-7  (Layne,  1975).  The  pitot  tube  is  constructed  of 
brass  and  is  mounted  in  the  third  window  of  the  test  combustor  0.64  cm 
from  the  disc  and  is  positioned  precisely  with  the  aid  of  a micrometer. 

The  static  pressure  is  measured  with  a pressure  tap  mounted  on  the  com- 
bustor wall  and  in  the  same  axial  plane  as  the  stagnation  pressure 
measurement.  The  static-stagnation  differential  pressure  is  monitored  in 
the  test  cell  with  a mercury  manometer  and  is  directly  related  to  the 
velocity  in  the  combustor  by  the  Bernoulli  equation  for  steady 
incompressible  flow. 

Since  the  pitot  tube  is  made  of  brass  and  is  not  water  cooled, 
only  cold  flow  measurements  were  possible.  The  results  of  the  measurements 
will  determine  if  these  cold  flow  data  are  extendable  to  velocities  during 
combustion.  Table  3-7  summarizes  the  operating  conditions  for  the  pitot 
tube  velocity  measurements.  These  measurements  were  also  made  with  the 
shrouded  combustor  under  the  first  condition  specified  in  Table  3-7  (see 
Figs.  4-1  and  4-2  for  the  results). 


Table  3-7.  Burner  Operating  Conditions  for 

Pitot  Ti4)e  Studies  (No  Combustion) 


m^  (kg/sec) 

p (atm) 

T.  (K) 
in  ^ ^ 

Vref  (m/sec) 

1.0 

2 

300 

25.4 

1.0 

4 

300 

12.7 

F.  Single  Point  Probe  Studies 

To  examine  the  prevaporizing  capability  of  the  fuel  preparation 
tube,  the  single  point  gas  sampling  probe  was  used  to  withdraw  san^les 
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from  flames  corresponding  to  three  different  combustor  operating 
conditions  using  two  different  nozzles.  These  operating  conditions  are 
listed  in  Table  3-8  and  represent  different  reference  velocities  and 
thus,  disparate  velocities  in  the  fuel  preparation  tube. 


Table  3-8.  Operating  Points  for  Single  Point  Probe  Studies 
(Fuel  = Jet  A,  ^ = 0.29) 


Nozzle 

Left 

Loac 

1/0 

jbvre  (1966) 
ling  Parameter 

f 

Air  Flow 

Rate 

m^  (kg/sec) 

d 

Inlet  Air 
Temperature 

T.  (K) 
in 

Combustor 

Pressure 

P (atm) 

^ ,,_1T75  2.75 
sec  atm  m 

16  and 
30-90 

.46 

1.0 

800 

8 

16  and 
30-90 

1.05 

1.0 

800 

5 

30-90 

1.55 

1.0 

800 

4 

16-90* 

2.57 

1.0 

800 

i 

3 

At  1.0  cm  downstream  of  the  disc  emission  concentration  data 
were  collected  with  the  probe  at  two  radial  positions  within  the  combustor 
(see  Fig.  3-8).  These  radial  positions  are  the  tube  centerline  and  the 
center  of  the  recirculation  zone  (i.e.,  halfway  between  the  inside 
diameter  of  the  tube  and  the  outside  diameter  of  the  disc) . Concentra- 
tions at  the  tube  centerline  are  inqjortant  for  examining  the  condition  of 
the  fuel  at  the  tube  exit  while  conditions  in  the  recirculation  zone 
provide  information  relevant  to  the  blowoff  analysis  which  is  a recircu- 
lation zone  phenomenon  (see  Tables  4-2  and  4-3).  Unfortunately,  combustor 
operating  conditions  proved  too  severe  for  single  point  tenperature 
measurements  in  these  regions. 


Lean  blowoff  occurred  at  this  condition  with  the  #30-90  nozzle. 


1 
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Probe  misalignment  errors  are  estimated  as  0.25  cm  (Tuttle 
et  al.,  1975)  in  radial  placement.  The  large  length  to  diameter  ratio  of 
the  probe  encourages  some  error  in  the  sample  tip  location  from  probe 
flexation  during  combustor  operation.  This  probe  tip  oscillation  about 
an  equilibrium  position  is  greatest  in  the  region  downstream  of  the  disc 
where  main  flow  streamlines  bend  toward  the  duct  centerline.  Hence,  in 
this  region  of  the  flow  field,  only  an  average  sanple  over  the  amplitude 
of  the  oscillation  is  possible. 

The  variety  of  experiments  mentioned  in  the  proceeding  sections 
are  fully  discussed  in  the  following  three  Chapters  and  provide  a guide 
to  the  development  of  the  characteristic  time  model  discussed  in 
Chapter  II. 
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IV.  FLOW  FIELD  CHARACTERIZATION 

Although  the  prevaporizing/premixing  geometry  discussed  in 
the  last  section  represents  an  extension  of  the  disc- in-duct  configu- 
ration studied  extensively  by  Tuttle  et  al.  (1976),  application  of  the 
characteristic  time  model  to  this  combustor  is  conplicated  by  the 
presence  of  the  fuel  preparation  tube  and  the  requirement  that  conplete 
fuel  vaporization  and  mixing  be  achieved  prior  to  combustion.  To  aid 
in  the  development  of  the  characteristic  time  model  for  combustion 
efficiency  and  flame  stabilization,  simple  flow  field  measurements  were 
made  in  the  region  of  the  fuel  preparation  tube.  These  measurements 
include  static  pressures  and  tenqieratures  along  the  tube  (inside  diameter) 
and  disc,  velocity  profiles  in  the  plane  of  the. disc  and  species  concen- 
trations at  the  tube  exit  and  in  the  recirculation  zone.  Particularly 
in^rtant  are  the  velocity  measurements  because  they  provide  information 
regarding  the  air  flow  split  between  the  tube  and  the  annulus  region  as 
well  as  the  average  equivalence  ratio  at  the  tube  exit.  Hopefully 
these  sinple  experiments  will  help  characterize  some  of  the  boundary 
conditions  in  the  flow  field. 

A.  Static  Pressure  and  Temperature 

For  the  instrumented  combustor  experiments,  static  pressures 
and  tenperatures  were  measured  in  the  regions  shown  in  Fig.  3-6  using 
the  shrouded  combustor  of  Fig.  3-4.  The  purpose  of  the  shroud  is  to 
smooth  the  incoming  flow,  shield  the  pressure  taps  and  thermocouples  and 
minimize  any  flow  disturbances  caiised  by  these  probes.  Table  4-1 
summarizes  the  data  obtained  in  cold  flow,  hot  flow  (preheat  combustor 
only)  and  hot  reacting  flow  (both  confcustors)  under  typical  combustor 
operating  conditions. 
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Overall  differential  pressure  measurements  (station  0-8)  show 
no  variation  between  cold  flow,  hot  flow  (preheat  combustor  only)  and 
reacting  flow  (both  combustors)  at  the  two  reference  velocities  tested. 
This  inplies  that  turbulence  levels  and  velocity  profiles  are  similar 
between  the  reacting  and  cold  flow  experiments  from  the  upstream  region 
to  the  disc.  Therefore,  cold  flow  velocity  measurements  in  the  plane 
of  the  disc  can  be  used  to  approximate  conditions  in  reacting  flow  at 
the  same  reference  velocity. 

Otiier  differential  pressure  measurements  are  very  low 
(<0.1  psid)  and  tend  to  strain  the  lower  range  of  the  pressure  trans- 
ducer. The  pressure  tap  located  in  the  ipstream  section  of  the  combustor 
(station  1)  measures  an  absolute  pressure  vdiich  is  used  as  reference  for 
the  other  readings;  the  remaining  three  are  differential  pressures  from 
this  reference  point.  One  result  is  evident;  the  pressure  at  station  4 
is  distinctly  higher  than  at  station  1 for  the  high  reference  velocity 
case.  This  indicates  that  some  degree  of  flow  separation  exists  at  the 
exit  of  the  fuel  preparation  tube.  Consequences  of  flow  separation  in 
the  mixing  tube  are  discussed  in  the  flashback  section  of  this  report. 

Thermocouples  on  the  inside  of  the  prevaporizing/premixing 
tube  (station  2 and  3)  provide  information  regarding  the  fuel  vaporiza- 
tion rate  at  the  wall.  For  the  reacting  flow  results,  the  teiiperatures 
(uncorrected)  measured  by  these  thermocoiples  are  substantially  lower 
than  the  combustor  inlet  tenperature.  This  temperature  decrease  may 
either  be  caused  by  the  liquid  Jet  A cooling  the  air  stream  as  it 
vaporizes  or  liquid  fuel  droplets  striking  the  thermocoiple  bead  and 
evaporating.  If  one  assumes  that  all  the  fuel  is  vaporized,  an  inlet 
tenperature  decrease  of  33K  is  anticipated  from  an  energy  balance. 
Consequently,  the  lower  tenperatures  reported  in  Table  4-1  are  probably 
a result  of  cooling  caused  by  droplet  evaporation  on  the  thermocouple 
bead.  For  the  low  velocity  experiments,  the  tenperature  at  station  3 
is  very  close  to  the  inlet  tenperature  indicating  better  vaporization 
near  the  wall  than  in  the  high  velocity  case. 


Finally,  the  disc  thermocouple  (station  6)  measures  the 
tenperature  of  the  gases  in  the  recirculation  zone  at  the  disc  surface. 
For  the  reacting  flow  experiments,  the  measured  temperature  is  approxi- 
mately 80?i  of  the  adiabatic  flame  temperature  based  on  the  overall 
equivalence  ratio.  It  is  the  recirculation  zone  that  provides  heat  and 
active  species  to  the  shear  layer  region  ^dlich  stabilizes  the  flame. 

The  purpose  of  the  hot  recirculation  zone  will  be  discussed  in  more 
detail  in  the  blowoff  section  of  Chapter  VI . 

B.  Velocity 

Results  of  the  overall  differential  pressure  measurements 
in  the  previous  section  show  that  velocities  at  the  tube  exit  for 
reacting  flow  can  be  estimated  by  performing  experimental  measurements 
under  cold  flow  conditions.  Utilizing  the  pitot  tube  of  Fig.  3-7  and 
the  operating  conditions  listed  in  Table  3-7,  velocities  in  a plane 
just  downstream  of  the  disc  were  obtained  at  two  different  reference 
velocities.  Profiles  are  shown  in  Fig.  4-1  and  Fig.  4-2  for  the  low 
and  high  reference  velocity  cases,  respectively,  using  both  the 
simplified  (Fig.  4-1  and  4-2)  and  shrouded  (Fig.  4-2)  combustor.  Note 
that  all  velocities  are  zero  near  solid  boundaries  except  at  the  housing 
wall.  This  is  due  to  the  combustor  windows  vdiich  are  recessed  in  the 
housing  and  are  in  the  same  plane  as  the  pitot  tube.  If  this  distance 
is  taken  into  account,  air  flow  rates  of  1.0  kg/ sec  are  predicted. 

In  both  cases,  the  velocity  at  the  tube  exit  can  be 
approximated  by  the  annular  velocity  V (velocity  based  on  the  area  of 
the  disc).  Experimentally  determined  air  flow  splits  are  shown  in  both 
figures  and  compare  favorably  with  the  calculated  value  of  m^  = 

0.275  kg/ sec  . The  velocity  profile  is  also  symmetric  about  the  com- 
bustor centerline;  in  addition,  nozzle  and  support  effects  are  dissipated 
in  the  fuel  preparation  tube  before  the  flow  reaches  the  tube  exit  as 
evidenced  by  the  flat  velocity  profiles  in  both  Fig.  4-1  and  4-2. 

The  major  effect  of  the  shrouded  conbustor  is  to  streamline 
the  flow  and  thus  flatten  the  velocity  profile  near  the  outer  region  of 
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x(an) 

Figure  4-2.  Cold  flow  pitot  tube  measurements  at  2 atm. 
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the  disc  (Fig.  4-2).  There  appears  to  be  only  a minimal  effect  on  the 
velocity  at  the  exit  of  the  prevapori zing/premixing  tube. 

C.  Species  Concentrations 

The  measurement  of  NO^  internally  in  a continuous  combustor 
provides  a useful  tool  for  the  verification  of  conputer  modeling  tech- 
niques. However,  in  this  combustor,  high  combustion  inefficiencies 
are  present  causing  large  amounts  of  unbumed  hydrocarbons  which  distort 
the  NO  emission  results:  this  effect  has  been  documented  in  Tuttle 

X ’ 

et  al.  (1974).  Also,  it  will  become  evident  later  in  this  chapter  that 
the  fuel  preparation  tube  has  failed  to  conpletely  prevaporize  the 
liquid  Jet  A.  For  these  reasons,  detailed  probing  of  the  combustor  was 
not  performed  in  this  preliminary  investigation;  rather,  single  point 
probe  measurements  were  made  at  the  tube  exit  (on  the  centerline)  and 
in  the  recirculation  zone  in  order  to  learn  more  about  fuel  prevapor- 
ization and  flame  stabilization,  respectively. 

Internal  measurements  were  obtained  under  the  burner  operating 
conditions  of  Table  3-8  at  two  radial  points  in  the  combustor  using  two 
different  nozzles.  Table  4-2  summarizes  the  results  for  measurements  on 
the  tube  centerline;  NO  emissions  are  not  tabulated  here  because  of 

X 

unavoidable  interactions  in  the  sampling  system.  Table  4-3  presents 

the  results  of  the  data  taken  on  the  centerline  of  the  recirculation 
zone  just  downstream  of  the  flameholder.  In  both  cases,  a high  velocity 
(low  P)  and  low  velocity  (high  P)  point  were  observed  along  with  the 
base  point. 

Examination  of  the  results  in  Table  4-2  show  that  there  is  no 
appreciable  nozzle  effect  on  either  the  unbumed  hydrocarbon  or  the  CO 
emissions  indexes.  Measured  HC  emissions  correspond  to  an  equivalence 
ratio  of  0.15-0.20  at  the  tube  exit.  This  is  very  close  to  the  overall 
equivalence  ratio  of  = 0.29;  however,  if  the  tube  equivalence  ratio 
is  confuted  based  on  the  air  flow  split  of  Section  B,  an  equivalence 
ratio  of  <))^  = 1.05  is  obtained  which  is  substantially  larger  than  the 
measured  value.  This  effect  is  probably  caused  by  the  failure  of  the 
fuel  preparation  tube  to  completely  vaporize  the  fuel.  If  liquid  drops 
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Table  4-2.  Results  of  Single  Point  Probe  Measurements 
on  Tube  Centerline  (0.64  cm  from  Disc) 

(m  . =1  kg/sec;  T.  = 800K;  tp  =0.29;  Fuel  = Jet  A;  Geometry  A) 

^ air  ^ in  ov  ’ 


Combustor  Pressure 

#16-90 

Nozzle 

#30-90  Nozzle 

P (atm) 

SMD  95ym 

SMD 

150uni 

HCEI 

COEI 

HCEI 

COEI 

8 

572 

14.8 

572 

21.6 

529 

17.0 

557 

17.0 

5 

629 

7.8 

572 

7.7 

572 

6.5 

629 

7.7 

4 

— 

— 

586 

10.0 

- - - 

— 

772 

9.3 

3* 


972  19.3 

858  14.8 


Table  4-3.  Results  of  Single  Point  Probe  Measurements 
in  Centered  Recirculation  Zone 
(0.64  cm  from  disc) 

(m  . =1  kg/sec;  T.  = 800K;  (p  =0.29;  Fuel  = Jet  A;  Geometry  A) 

air  * in  ov 


Combustor  Pressure 
P (atm) 

#16-90  Nozzle 
95ym 

NO  El  HCEI 

X 

COEI 

#30-90  Nozzle 
SMD  ''j  150ym 
NO  El  HCEI 

X 

COEI 

8 

5.4 

0.2 

1.8 

6.8 

2.9 

4.0 

6.4 

1.5 

2.4 

6.2 

0.4 

4.0 

3.5 

0.2 

6.1 

4.3 

1.1 

8.6 

4 4 

4.5 

4.1 

0.5 

12.4 

3.B 

1.3 

16.1 

IS.' 

^ » )« 
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are  present  at  the  tube  exit,  they  (depending  on  their  size)  may  not 
be  drawn  into  the  sanpling  probe  and  measured  by  the  hydrocarbon 
detector.  A non-uniform  fuel  distribution  can  also  be  responsible 
for  the  relatively  low  equivalence  ratios  observed  on  the  tube  center- 
line.  High  fuel  injection  pressures  and  large  droplet  evaporation 
times  may  force  large  quantities  of  fuel  to  the  tube  wall  and  not 
allow  sufficient  time  for  prevaporization  and  premixing  prior  to 
combustion. 

Recirculation  zone  conditions  (Table  4-3)  show  low  unbumed 
hydrocarbon  and  carbon  monoxide  with  relatively  high  NO^  emission.  This 
inplies  high  temperatures  (see  Section  A)  and  long  mixing  times  favorable 
to  NO  formation  and  CX)  and  HC  oxidation.  Again,  no  appreciable  fuel 
injector  effect  is  observed  except  for  a slight  increase  in  CO  emissions 
with  increasing  droplet  34D. 


D.  Summary 

It  has  been  demonstrated  in  Section  A of  this  chapter  that  for 
high  reference  velocity  flows,  a distinct  region  of  flow  separation 
exists  within  and  near  the  end  of  the  prevaporizing  tube.  Also,  there 
is  a great  possibility  that  fuel  droplets  are  reaching  the  tube  exit  (at 
least  near  the  tube  wall)  as  indicated  by  static  temperature  measurements. 

However,  equally  important  is  the  fact  that  in  all  cases  tested 
the  overall  differential  pressure  drop  between  the  combustor  inlet  and 
the  disc  is  the  same  for  cold  flow,  hot  flow  (preheat  combustor  only) , 
and  reacting  flow  (both  combustors) . This  means  that  velocity  measure- 
ments at  the  disc  can  be  made  under  cold  flow  conditions  and  related  to 
a reacting  flow  environment.  Using  these  results,  it  was  found  that  the 
air  velocity  in  the  fuel  preparation  tube  can  be  estimated  from  the 
velocity  confHJted  on  the  basis  of  the  disc  blockage.  This  result  will 
prove  valiiohlc  in  the  next  chapters  for  estimating  not  only  velocities 
’t  ■.  ratio*  at  the  tii>e  exit. 

~ • ' t ■ f i * 


t tt>» 


70 


however,  HC  concentrations  are  not  as  high  as  one  would  expect  if  the 
fuel  were  completely  vaporized  and  mixed  with  the  air.  In  the  recircu- 
lation zone,  the  opposite  is  true;  low  CX)  and  HC  emissions  as  well  as  high 
temperatures  are  present.  This  is  in  agreement  with  Zukoski  and  Marble 
(1956)  who  state  that  the  recirculation  zone  is  a region  of  hot  burned 
gases  providing  sufficient  thermal  energy  for  ignition. 

The  following  chapters  are  concerned  primarily  with  input/output 
experiments  rather  than  detailed  single  point  measurements.  Both  com- 
bustion efficiency  and  flame  stabilization  are  considered  using  the 
geometry  of  Fig.  3-3.  In  addition,  the  degree  of  prevaporization  is 
varied  by  injecting  two  different  fuels.  Jet  A and  liquid  propane. 
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V.  OCMBUSTiraj  EFFICIENCY 


A.  Quantification  of  Combustion  Efficiency 
Combustion  efficiency,  n^.,  is  defined  as  the  percentage  of 
available  chemical  energy  released  i5X)n  combiistion  of  a fuel  and  oxidizer. 
One  hundred  percent  combustion  efficiency,  however,  is  not  attainable 
due  to  equilibrium  conditions  (especially  at  elevated  temperatures) ; also 
inccOTplete  combustion  through  quenching  can  limit  combustion  efficiency. 

This  unused  chemical  energy  is  a measure  of  combustion  inefficiency, 
and  manifests  itself  primarily  in  concentrations  of  carbon  monoxide  (CO) 
and  unbumed  hydrocarbons  (HC) . Hence  the  Air  Force  equation  (Blazowski 
and  Henderson,  1974)  readily  quantifies  comb'ostion  inefficiency  from  practi- 
cal systems  through 


1-n, 


(EICO)  + (EIHC) 


'‘fuel 


10^ 


5-1 


where  EIX  = emissions  index  of  species  X (g  X/kg  fuel) 
Q = respective  lower  heating  value  (cal/g) 
and  HC  are  measured  as  in  the  present  study. 


B.  Combustion  Efficiency  Correlation 
Several  parameters  are  available  which  correlate  combustion 
efficiency  from  reacting  systems.  These  include  percent  of  maximum 
rated  power  (Hammond,  1977),  combustor  reference  velocity  (Grobman,  1972), 
and  air  loading  (Longwell  and  Weiss,  1955).  A traditional  correlation 
for  combustion  efficiency  in  gas  turbine  combustors  is  the  9 parameter 
(Lefebvre,  1966).  This  parameter  is  defined  in  Eq.  3-1  but  reiterated 
here  for  convenience 

P • coiriNiiitor  prv««ur«,  atm 
T._  • inlat  atr  ta^rraturv,  t 

I * 


5-2 
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m„  = air  mass  flow  rate,  kg/sec 

® 2 
A r = maximum  combustor  cross  sectional  area,  m 
ref 

= maximum  ccanbustor  width,  m. 

The  6 parameter  is  a powerful  design  tool  in  that  for  a given 
combustor,  it  will  collapse  efficiency  data  onto  a single  curve  for  any 
combination  of  airflow,  inlet  tenperature  and  pressure.  Unfortunately, 
the  correlation  does  not  account  for  explicit  variations  in  gecsnetry 
(for  instance  a change  in  the  axial  location  of  the  secondary  air  pene- 
tration jets  in  a conventional  burner),  equivalence  ratio,  fuel  injector 
or  fuel  type  (Dodds  et  al.,  1975).  Correlating  deficiencies  in  the  6 
parameter  do  however  help  identify  relative  trends,  particularly  for 
a combustor  exhibiting  a radical  departure  in  design  from  conventional 
configurations  (such  as  the  tube-and-disc  burner  of  this  study) . Further- 
more, the  6 parameter  can  demonstrate  when  these  effects  are  insignificant. 

Combustion  efficiency  data  from  the  baseline  (Geometry  A)  tube- 
and-disc  flameholder  are  illustrated  as  a function  of  6 in  Fig.  5-1. 

Here  data  for  two  simplex  pressure  atomizing  nozzle  (SMD)  variations 
utilized  in  the  Jet  A experiments  are  presented  at  a constant  inlet  air 
temperature  of  BOOK,  equivalence  ratio  of  0.29  and  air  flow  rate  of  1.0 
kg/s  (recall  from  Table  3-3  that  these  and  a combustor  inlet  pressure 
of  5.0  atm  define  the  base  operating  point).  It  appears  from  this  figure 
that  heterogeneous  effects  on  combustion  efficiency  are  negligible  for 
this  geometry,  these  flow  conditions,  and  this  nozzle  range  except  perhaps 
at  low  0 where  nozzle  effects  are  typically  more  pronounced.  Droplet 
effects  upon  efficiency  should  become  more  important  at  the  lower  inlet 
air  temperatures,  lower  equivalence  ratios,  and  lower  combustor  pressures 
(Lefebvre,  1966). 

These  results  also  illustrate  the  very  low  combustion  efficiencies 
associated  with  this  type  of  burner  especially  at  high  reference  velocities 
(low  0).  Dodds  et  al.  (1975)  measured  combustion  efficiencies  from  the 
disc- in-duct  burner  discussed  in  Chapter  II  in  the  same  manner  and  found 
■ininun  efficiencies  of  85  percent  at  a 6 value  of  0,10.  Thus  the  location 
of  the  fuel  injector  with  respect  to  the  recirculation  zone  has  a dramatic 
e<fe»»  '^wistion  effKiencv,  a tipiificant  finding  relevant  to  the 
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the  tube-and-disc  configuration  for  both  Jet  A and  propane  (see  also  Fig. 

5-3). 

The  effect  of  tube-and-disc  flameholder  configuration  is  illustrated 
in  Fig.  5-2  for  Jet  A and  the  constant  base  point  conditions  enumerated 
previously.  Geometry  effects  are  most  important  in  the  low  to  mid  0 regime 
where  in  order  of  decreasing  combustion  efficiency  the  geometries  order 
as  D or  E (small  diameter  fuel  preparation  tube),  B,  A and  C (large  diameter 
tube) . High  6 or  high  pressure  efficiencies  progress  toward  geometry 
independence  due  to  more  effective  kinetics  and  reduced  dissociation  losses, 
better  fuel  spray  evaporation,  and  increased  residence  (or  mixing)  times 
within  the  combustor.  Only  three  of  the  five  flameholder  geometry  varia- 
tions were  tested  with  propane.  Combustion  efficiencies  from  these  con- 
figurations are  displayed  in  Fig.  5-3  at  the  identical  operating  conditions 
of  Fig.  5-2.  These  three  configurations  exhibit  the  same  relative  effi- 
ciencies as  demonstrated  with  Jet  A.  That  is,  the  small  tube  flameholder 
is  more  efficient  than  the  large  tube  variations  with  geometry  B again 
exhibiting  higher  efficiencies  than  A. 

The  effect  of  inlet  air  tenperature  at  a constant  equivalence 
ratio  and  combustor  pressure  of  0.29  and  5 atm  is  displayed  in  Fig.  5-4 
for  the  base  geometry  burning  propane  and  Jet  A.  As  expected  an  increase 
in  inlet  air  temperature  enhances  the  rate  of  fuel  oxidation  to  CO  and  of 
CO  to  CO2  thus  increasing  combustion  efficiency. 

Obvious  from  Fig.  5-4  and  illustrated  in  Fig.  5-5  is  the  effect 
of  fuel  type  on  combustion  efficiency  in  the  baseline  burner.  Jet  A shows 
consistently  higher  combiistion  efficiencies  than  propane,  a result  imexpected 
from  a viscosity,  volatility,  heating  value  or  fuel  state  argument  (see 
Lefebvre,  1966  or  Tuttle  et  al.,  1977).  Note  that  efforts  to  include  fuel 
type  in  the  efficiency  correlation  are  precluded  since  Jet  A has  a lower 
heat  of  combustion  relative  to  propane;  that  is  Eq.  5-1  predicts  higher 
efficiencies  from  propane  rather  than  Jet  A. 

Equivalence  ratio  also  plays  a significant  role  in  determining 
combustion  efficiency  in  the  tube-and-disc  flameholder.  Specifically, 

Fig.  5-6  shows  that  for  both  propane  and  Jet  A increasing  the  overall 
equivalence  ratio  at  constant  combustor  loading  causes  combastion  efficiency 
to  deteriorate.  This  trend,  whi  h (vcurreil  for  all  gismietrv  and  furl  varia 


Figure  5-2.  Tube -and -disc  flameholder  geometry  effect  on 
combustion  efficiency  (Jet  A) . 


Figure  5-3.  Tube-and-disc  flameholder  geometry 
effect  on  combustion  efficiency 
(Propane) . 


Figure  5-6.  The  overall  equivalence  ratio  effect 
on  combustion  efficiency  in  the  tube- 
and-disc  burner. 
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tions,  is  opposite  that  observed  in  the  disc-in-duct  (Tuttle  et  al.,  1977) 
and  conventional  (Mellor,  1977b)  combustors  for  essentially  the  same  range 
in  overall  equivalence  ratio  (i.e.,  0.1  0.4).  If  this  effect  were 

observed  for  Jet  A only  (i.e.,  a negligible  or  opposite  equivalence  ratio 
trend  on  combustion  efficiency  for  propane)  droplet  effects  would  be  siis- 
pect;  but  such  is  not  the  case  here. 

Fig.  5-5  in  conjunction  with  Fig.  5-6  yields  valuable  insight 
into  the  phenomena  which  appear  to  affect  combustion  efficiency  in  the 
tube-and-disc  combustor.  First,  the  eqiri valence  ratio  effect  where  de- 
creasing efficiencies  acconpany  increases  in  overall  equivalence  ratio 
suggests  that  equivalence  ratios  evaluated  in  the  preparation  tube  are 
important  for  combustion  efficiency.  Specifically,  most  tube-and-disc 
geometries  produce  fuel  rich  equivalence  ratios  in  the  tube  throughout 
the  combustor  operating  range  (except  configurations  A and  B which  have 
one  lean  tube  equivalence  ratio  of  0.85).  Thus,  increasing  the  equivalence 
ratio  in  the  tube  (<ti  > 1)  lowers  the  burned  gas  temperature  and  the 
combiistion  efficiency.  The  importance  of  the  tube  equivalence  ratio  to 
combustion  efficiency  is  not  unreasonable  upon  consideration  of  a true 
dunp  configuration  where  air  is  only  supplied  to  the  combustor  through 
the  fuel  preparation  tube.  In  addition,  rich  equivalence  ratios  in  the 
tube  may  explain  why  higher  efficiencies  are  obtained  with  Jet  A than 
with  propane.  Since  propane  flash  vaporizes  upon  injection,  the  flame 
is  siq)plied  with  a mixture  having  a richer  vapor  phase  equivalence  ratio 
than  Jet  A;  thus,  higher  combustion  efficiencies  are  obtained  with  Jet  A 


^'^tube 


closer  to  stoichiometric) . 

Problems  related  to  the  sampling  system  may  also  explain  the 


observed  fuel  type  effects  on  combustion  efficiency.  First,  it  is  likely 


that  the  distillate  nature  of  Jet  A promoted  some  HC  condensation  within 


the  sanpling  system.  Dodds  et  al.  (1975)  attributed  the  higher  combustion 


efficiencies  realized  from  #5  fuel  oil  than  from  Jet  A to  this  condensation 


phenomenon.  Secondly,  since  the  sanpling  systan  is  designed  to  induct 
gaseous  emissions,  any  liquid  Jet  A droplets  present  at  the  exit  plane  will 
probably  bypass  the  sanpling  rake.  Both  of  these  effects  would  promote 
deceptively  higher  combustion  efficiencies  (lowr  HC  readings)  with  Jet 
A relative  to  propane. 
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The  present  effort  in  correlating  combustion  efficiency  in 
the  tube -and -disc  flameholder  is  devoted  to  extending  the  characteristic 
time  model  (discussed  in  Chapter  II)  to  this  configuration.  As  noted 
previously,  inclusion  of  fuel  type  in  these  correlations  may  require 
more  probing  experiments  based  on  the  relative  combustion  efficiencies 
of  Jet  A and  propane.  However,  it  is  expected  that  the  characteristic 
time  model  will  correlate  inlet  air,  equivalence  ratio  and  geometry 
effects  for  at  least  the  propane  data.  In  particular.  Fig.  5-4  and  5-5 
indicate  that  heterogeneous  effects  are  not  significant  to  combustion 
efficiency  in  the  tube-and-disc  burner  within  the  range  of  variables 
incorporated.  Therefore,  the  CO  correlation  of  Tuttle  et  al.  (1977)  as 
expressed  in  Eq.  2-24  should  extend  to  combustion  inefficiency  (provided 
that  HC  are  proportional  to  CO)  giving 

100-n^  'V  5-3 

where  = chemical  kinetic  time  for  fuel  and  CO  oxidation 
to  CO2  (Arrhenius  form) 

Xg^  = turbulent  mixing  time  in  the  region  of  predominant 
CO  and  HC  quenching  (generally  a macro  length  scale 
divided  by  a velocity) . 

Inlet  air  temperature  effects  will  be  included  in  the  characteristic 
time  correlation  primarily  through  the  chemical  kinetic  time  for  fuel  oxidation, 
x^.  Inclusion  of  equivalence  ratio  is  also  expected  through  this  kinetic 
time  by  utilizing  the  burned  gas  temperature  based  on  the  total  tube 
equivalence  ratio  and  by  accounting  for  inlet  air  vitiation.  The  turbulent 
mixing  time,  x^j^,  should  permit  inclusion  of  flameholder  geometry  effects 
through  proper  determination  of  the  length  and  velocity  scales.  However, 
the  appropriate  characteristic  times  for  combustion  efficiency  have  not 
yet  been  determined  quantitatively. 

Combustion  inefficiency  was  also  measured  with  the  shrouded 

combustor  of  Fig.  3-4.  Results  are  sumarized  in  Table  5-1  for  the  burner 

operating  conditions  of  Table  3-5.  In  general,  the  shrouded  con4>ustor  is 

slightly  more  efficient  than  the  sinfilified  coni^ustor.  Since  the  major 

parameters  which  influence  efficiem-v  fV’  r.  T.  , etc.),  hav-e  not  been 

reT  in 

chanited  these  results  *re  tw>t  taiexpes  te\l . Therefore,  ^isnuiter  modeling 
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Table  5-1.  Effect  of  Slirouded  Combustor  on  Combustion 
Inefficiency  ((J)ov=0.29,  Jet  A) 


Combustor 

Pressure 

P(atm) 

Inlet 

Tarperature 
Tin  (K) 

Reference 

Velocity 

Vref  (m/sec) 

Ccanbustion 

Inefficiency 

lOO-fiQ 

Combustor 

8 

800 

17.0 

.02-. 03 
.02-. 06 

Sinplified 

Shrouded 

5 

800 

27.1 

.17-. 18 
.14-. 15 

Simplified 

Shrouded 

4 

800 

34.0 

.42-. 43 
.37-. 44 

Simplified 

Shrouded 

for  the  shrouded  combustor  should  demonstrate  only  a second  order  effect 
of  the  shroud. 

In  order  to  obtain  better  understanding  of  the  flovrfield,  de- 
tailed internal  measurements  are  proposed  for  selected  tube-and-disc 
geometries.  Single  point  probe  measurements  will  provide  species  con- 
centration profiles  and  temperature  estimates  and  should  yield  insight 
into  the  most  significant  region(s)  of  combustion  inefficiency,  be  it  in 
the  shear  layer  encompassing  the  recirculation  zones  or  some  region  down- 
stream. In  addition,  velocity  magnitude  and  direction  measurements  via 
laser  Doppler  velocimetry  are  expected  to  further  aid  in  the  identification 
of  appropriate  characteristic  time  parameters. 


83 


h 

I 


VI.  FLAME  STABILIZATION 


One  obvious  requirement  of  advanced  prevaporizing/premixing 
combustors  is  proper  flame  stabilization.  High  approach  flow  velocities 
and  relatively  low  fuel  flow  rates  may  cause  the  flame  to  be  extinguished 
(blowoff)  under  certain  low  power  operating  conditions.  This  is 
particularly  inportant  in  prevaporizing/premixing  combustors  because 
locally  high  equivalence  ratios  and  temperatures  normally  present  in 
conventional  combustors  are  absent.  Thus,  lean  blowoff  as  well  as 
ignition  limits  are  much  narrower. 

Introduction  of  the  fuel  upstream  of  the  flameholder  also 
creates  a second  flame  stabilization  problem  known  as  flashback,  which 
can  cause  considerable  damage  to  the  combustion  chamber.  Flashback 
results  from  increasing  the  fuel/air  ratio  toward  stoichiometric 
(increasing  the  flame  speed)  or  decreasing  the  air  velocity,  allowing 
the  flame  to  propagate  into  the  fuel  preparation  tube. 

In  the  following  sections,  problems  associated  with  flashback 
and  blowoff  are  examined  for  both  laminar  and  turbulent  flow  situations. 
With  the  help  of  existing  literature,  the  characteristic  time  model 
discussed  in  Qiapter  II  is  successfully  applied  to  lean  blowoff  in 
turbulent  bluff  body  stabilized  flames.  In  addition,  a qualitative 
model  has  been  developed  to  explain  the  reported  upstream  flame 
propagation.  Both  these  results  are  extended  to  flame  stabilization 
data  obtained  using  the  simplified  pre vaporizing/premixing  combustor 
described  in  Chapter  III.  The  blowoff  model  for  the  tube-and-disc 
combustor  correlates  variations  in  pressure,  inlet  temperature, 
reference  velocity,  flameholder  size  and  degree  of  prevaporization. 


A.  Derivation  of  the  Correlation 
In  the  previous  sections,  the  characteristic  time  model  was 
applied  to  both  practical  burners  and  combustors  of  research  interest 
for  the  purpose  of  emission  control.  However,  flame  stabilization  can 
also  be  interpreted  via  the  characteristic  time  approach.  For  laminar 
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premixed  gaseous  (Bunsen  burner  type)  flames,  the  model  is  essentially 
equivalent  to  the  Peclet  number  correlation  proposed  by  Putnam  and  Jensen 
(1949) ; it  also  represents  a quantitative  application  of  the  Zukoski  and 
Marble  (1956)  theoiy  for  predicting  lean  blowoff  limits  in  turbulent 
bluff  body  stabilized  flames.  In  all  cases,  flame  stabilization  involves 
the  conqjetition  between  fluid  mechanic  and  chemical  effects . The  model 
therefore  characterizes  the  combustion  process  by  an  appropriate  fluid 
mechanic  and  chemical  time  scale  evaluated  in  the  particular  I'egion  of 
the  flow  field  of  importance  to  the  stabilization  process. 

1.  Laminar  Premixed  Flames 

The  majority  of  flashback  and  blowoff  data  available  in  the 
literature  has  been  obtained  in  laminar  premixed  flows  for  gas -burning 
appliance  applications  at  room  temperature  and  atmospheric  or  sub- 
atmospheric  pressure  (Putnam  and  Jensen,  1949;  Wohl  et  al.,  1949;  Grumer 
et  al.,  1953  and  1956;  Khitrin  et  al.,  1965).  Turbulent  approach  flows 
(Wohl  et  al.,  1949;  Wohl,  1953;  Grumer  et  al.,  1953  and  1956;  Boyer  and 
Friebertshauser,  1957;  Fine,  1958;  Grumer,  1958;  Khitrin  et  al.,  1965) 
and  geometry  variations  (Grumer  et  al.,  1953  and  1956)  have  also 
been  considered  under  these  conditions.  Only  a few  investigators 
(Dugger,  1950;  Grumer  et  al.,  1956;  Khitrin  et  al.,  1965)  have  reported 
flashback  and  blowoff  results  for  approach  flows  at  elevated  tenperatures 
(see  Table  6-1);  however,  all  are  for  gaseous  premixed  systems  burning 
H2  or  pure  hydrocarbons.  After  a brief  review  of  existing  correlations, 
some  of  the  data  listed  in  Table  6-1  will  be  used  to  illustrate  applica- 
tions of  the  characteristic  time  model  and  also  to  demonstrate  the 
equivalence  of  the  model  to  the  Peclet  nunter  correlation. 

For  laminar  flow  in  a tube,  flame  stabilization  involves  the 
conpetition  between  fluid  mechanics  and  chemical  kinetics  in  the  boundary 
layer.  Early  work  on  flame  stabilization  in  laminar  flow  consisted  of 
experimentally  determining  the  average  velocity  at  which  flashback  or 
blowoff  occurred  as  a function  of  equivalence  ratio.  However  results  of 
this  nature  were  only  valuable  for  one  particular  geometry  and  inlet 
condition. 
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Table  6-1.  Flashback  and  Blowoff  Data  for 

Heated  Reactants  (Laminar  Flow,  P*1  atm) 


Limit 

Fuel/Oxidizer 

Inlet 

Temperature,  K 

Equivalence 

Ratio 

Reference 

Blowoff 

CjHj/Air 

300  - 620 

0.8  - 

1.6 

Dugger  (1950) 

Flashback 

CH^/Alr 

300  - 525 

0.6  - 

1.4 

Grumer  et  al. 
(1956) 

Blowoff 

CH./Alr 

300  - 525 

0.6  - 

1.4 

Grumer  et  al. 
(1956) 

Flashback 

CH/Oj 

373  - 673 

0.3  - 

1.6 

Khitrin  et  al 
(1965) 

El 


The  velocity  gradient  theory  of  Lewis  and  von  Elbe  (see 
disciission  in  Lewis  and  von  Elbe,  1961)  represents  one  of  the  first 
attempts  at  correlating  flame  stabilization  data.  For  flashback,  the 
theory  balances  the  velocity  gradient  at  the  wall  (g)  with  the  ratio 
of  laminar  flame  speed  (S, ) divided  by  the  quenching  distance  (d  ) , 

Lt  C 


= ^ 
®fb  dr 


wall 


6-1 


assuming  that  the  stream  velocity  curve  can  be  replaced  by  a straight  line 
having  the  slope  of  the  stream  velocity  curve  at  the  wall.  If  the 
velocity  gradient  is  less  than  S, /d  , the  flame  will  flash  back;  if  it 
is  greater  than  S.  /d  , the  flame  will  blov  out  of  the  tube  and  stabilize 
on  the  burner  rim.  Blowoff  is  represented  in  a similar  maimer  except 
the  velocity  gradient  at  the  wall  is  replaced  by  the  gradient  at  the 
boundary  between  the  tube  flow  and  outside  ambient  conditions.  When 
gj^  exceeds  the  critical  value  of  Sj^/d^  at  the  free  boundary,  the 
flame  blows  off.  Therefore  plots  of  g versus  equivalence  ratio  (4)) 
should  indicate  regions  of  stable  and  unstable  combustion.  A typical 
diagram  of  velocity  gradient  as  a function  of  equivalence  ratio  is 
shown  in  Fig.  6-1.  Note  that  in  air  the  blowoff  curve  continues  to  rise 
after  the  stoichiometric  equivalence  ratio  is  reached;  this  is  due  to  the 
diffusion  of  oxygen  into  the  rich  mixture  wfiich  increases  the  local  flame 
speed. 


fuel  lean  <p=l  fuel  rich 


Figure  6-1.  Representation  of  flashback  and  blowoff 
limits  for  a laminar  premixed  hydrocarbon 
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The  velocity  gradient  at  the  wall  is  easily  catputed  in  laminar 
flow  by  using  the  friction  factor  (f)  defined  by  Eq.  6-2, 


wall 


u'^  = 8vg  vT^ 


6-2 


where  v = kinematic  viscosity 
and  u = mean  velocity  of  the  approach  flow. 
By  combining  Eq.  6-1  and  6-2,  the  velocity  gradient  becomes 


g = 


a-  - 8v  • 
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Although  the  velocity  gradient  correlation  (g  vs.  4))  includes 
variations  in  tube  diameter,  it  fails  to  account  for  changes  in  pressure, 
temperature  and  inlet  conposition.  In  1949,  Putnam  and  Jensen  proposed 
the  Peclet  number  correlation  for  flashback  which  results  in  a more 
universal  curve  rather  than  a family  of  curves  dependent  upon  inlet 
conditions.  Assuming  the  quenching  distance  to  be  proportional  to 
the  thermal  diffusivity  (a)  divided  by  the  flame  speed,  Eq.  6-3  becomes 


6-4 


For  laminar  flow,  the  friction  factor  is 

64 


f = 


Re 
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Substituting  Eq.  6-5  into  Eq.  6-4  and  multiplying  by  d^/a  (d  = tube 
diameter)  gives 


a a 


6-6 


or 


Pec  '''  Pe^ 
Sl  u 


6-7 
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where  Peg  = Peclet  tiumber  based  on  flame  speed 

and  Pe^^  * Peclet  number  based  on  u. 

Putnam  and  Jensen  (1949)  and  Khitrin  et  al.  (1965)  have  been  able  to 
successfully  correlate  laminar  flasliback  in  tubes  using  Eq.  6-7.  An 
exanple  of  the  Peclet  nunfcer  correlation  is  shown  in  Fig.  6-2  for  the 
flashback  data  of  Khitrin  et  al.  (1965)  at  elevated  inlet  ten?)eratures 
(see  Table  6-1);  common  symbols  represent  data  points  corresponding  to 
inlet  tenperature  Arariations  at  a given  equivalence  ratio. 

Although  the  Peclet  number  correlation  is  more  general  than 
plots  of  g versus  0,  it  requires  accurate  determination  of  the  laminar 

flame  speed  vdiich  is  not  always  possible.  An  alternate  approach  is  to 

view  the  flame  stabilization  process  in  terms  of  characteristic  times. 
Since  the  reaction  rate  is  related  to  the  laminar  flame  speed 

(reaction  rate)*^  , 6-8 

an  ignition  delay  time  can  be  defined  as 

-1  -2 

'V/  (reaction  rate)  'v  . 6-9 


Furthermore,  assuming  that  the  velocity  gradient  at  the  wall  (flashback) 
or  free  boundary  (blowoff)  is  proportional  to  a fluid  mechanic  time 
for  laminar  flow,  Eq.  6-10  becomes 


6-10 


Thus,  Eq.  6-4  involves  the  competition  between  a fluid  mechanic  (t^j^) 
and  chemical  (tj^^)  time 

^fm  ^hc 

at  the  limit. 

The  fuel  ignition  delay  time  is  computed  from  the  general 
overall  rate  expression  for  hydrocarbon  oxidation 


Figure  6-2.  Peclet  number  correlation  for  014/02 

flashback  data  (Khitrin  et  al.,  1965) 

at  various  (f)  and  T.  . 

in 
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^ ■V  [HCl  [Ojl  (2^) 


6-12 


assuming  a first  order  dependence  of  [HC]  and  [O2]  (overall  second  order) . 
Since  most  practical  flows  of  interest  are  fuel  lean,  the  [O2]  concentration 
should  remain  approximately  constant  so  that  Eq.  6-12  reduces  to 
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Dividing  by  the  overall  density  of  the  system,  Eq.  6-13  in  terms  of  mole 
fractions  becomes 
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which  has  units  of  sec  Defining  the  ignition  delay  time  as 
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yields 
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since  the  approach  flow  equivalence  ratio  is  directly  proportional  to  the 
hydrocarbon  mole  fraction. 

Combining  Eq.  6-16  and  Eq.  6-10  with  Eq.  6-11,  the  characteristic 
time  model  for  flame  stabilization  in  laminar  flow  becomes 
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Note  that  Eq.  6-17  is  a semi-empirical  expression  since  the  activation 
energy  (E)  is  chosen  to  give  the  best  fit  to  the  experimental  data.  The 
flashback  data  reported  by  Khitrin  et  al.  (1965)  are  now  replotted  in 
Fig.  6-3  using  this  equation.  The  appropriate  hydrocarbon  ignition  delay 
time  was  found  to  be 


^ ^q-7  exp  (18000/RT) 
he  ♦ 


6-18 


1 2 5 10  20  50  100  200 


X 10^  (sec) 


Figure  6-3.  Characteristic  time  model  for  CHa/Oo 
flashback  data  (Khitrin  et  al.,1965j 
at  various  and  T.  . 
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(units  on  E are  cal/nxDle)  with  T evaluated  at  the  average  between  the 
inlet  and  the  adiabatic  flame  tenperature  at  flashback.  However,  the 
apparent  activation  energy  in  Eq.  6-18  is  substantially  lower  than 
values  reported  in  Cooke  and  Williams  (1975)  for  methane.  This  is 
attributed  to  the  selection  of  an  average  temperature  and  first  order 
fuel  dependence  in  Tj^^.  The  characteristic  time  model  of  Fig.  6-3  is 
extremely  accurate  for  equivalence  ratios  less  than  unity  (rich  flash- 
back data  are  not  expected  to  correlate  based  on  the  derivation  of 
with  a correlation  coefficient  very  close  to  one.  Therefore,  the  Peclet 
number  and  characteristic  time  correlations  for  flashback  in  laminar 
flow  are  essentially  identical  as  the  derivation  predicts,  for  equivalence 
ratios  of  practical  interest. 

At  the  present  time,  very  few  flashback  data  are  available  for 
practical  combustors.  Thus,  tremendous  extrapolation  of  available  flame 
stabilization  data  in  terms  of  fuel,  pressure,  inlet  tenperature,  geometry 
and  velocity  typical  of  continuous  combustors  is  required.  Altenkirch 
and  Mellor  (1973)  found  it  difficult  to  compare  their  lean  blowoff  data 
with  Spalding's  (1955)  Peclet  number  correlations  because  the  laminar 
flame  speed  in  PcSl  unknown  even  with  C^Hg  at  4 atm  and  950K  inlet 
tenperature . For  flashback  there  is  a question  as  to  vdiether  it  is 
actually  a boundary  layer  phenomenon,  especially  in  turbulent  flows. 
Flashback  has  been  observed  to  occur  on  the  flow  axis  in  nozzles  with 
relatively  flat  velocity  profiles  (Wohl,  1953)  and  for  turbulent  flows 
in  large  diameter  tubes  (Grumer,  1958).  Furthermore,  heterogeneous 
effects  such  as  droplets  may  also  be  inportant.  Because  of  these 
difficulties,  the  authors  believe  the  characteristic  time  model  provides 
a greater  degree  of  flexibility  for  incorporating  these  types  of  phenomena 
in  flame  stabilization  correlations  for  modem  air-breathing  engine 
applications. 

2.  Prevapori zing/ Premixing  Combustors 

Since  laminar  flame  speeds  required  for  both  the  velocity 
gradient  and  Peclet  number  correlation  are  extremely  difficult  to 
calculate  for  conditions  typical  of  modem  aircraft  combustors,  we  prefer 
a coiTelation  technique  vdiich  is  directly  related  to  inlet  conditions. 

The  inverse  theta  parameter  (9  ^)  of  Lefeb'vre  (1966)  is  comnonly  used 
to  correlate  flame  stabilization  data  (and  combustion  efficiency  - see 
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Chapter  V)  in  gas  turbines.  6 which  represents  combustor  loading, 
is  given  by 
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where  m^  is  the  air  flow  rate,  P combustor  pressure,  inlet  air 
temperature  and  combustor  area  (diameter).  Fig.  6-4  illus- 

trates the  relationship  which  exists  between  the  classical  flashback  and 
blowoff  limits  for  a fuel  lean  combustion  system.  Flashback  is  most 
prevalent  under  conditions  of  high  local  equivalence  ratio  and  relatively 
low  combustor  loading  (low  reference  velocity) . Note  that  for  large 
0 ^ a point  is  reached  where  flashback  is  not  possible  over  any  range 
of  local  equivalence  ratios.  Since  blowoff  is  meclianistically  similar 
to  flashback,  the  lean  blowoff  limit  is  characterized  by  relatively 
high  velocities  and  low  equivalence  ratios. 

In  the  following,  flashback  and  blowoff  in  prevaporizing/ 
premixing  combustors  are  discussed.  A recent  review  (Plee  and  Mellor, 
1977)  of  "flashback"  reported  in  the  literature  has  identified  four 
mechanisms  which  may  lead  to  upstream  flame  propagation.  In  addition, 
a characteristic  time  model  has  been  developed  for  correlating  lean 
blowoff  data  from  bluff  body  stabilized  flames  burning  gaseous  propane. 


a.  Flashback 

For  realistic  advanced  prevaporizing/premixing  combustors,  no 
con^irehensive  study  involving  flashback  is  currently  available  in  the 
literature  for  practical  hydrocarbon  fuels  at  elevated  inlet  temperatures 
and  pressures.  Instead,  iqistream  flame  stabilization  is  usioally 
observed  while  performing  combustion  efficiency  or  emission  measurements 
or  during  transient  combustor  operation.  A great  deal  of  overlap  and 
inconsistency  exists  between  these  reported  flashbacks  and  the  stable 
operating  conditions  of  other  investigators  vdiich  cannot  be  interpreted 
by  single  inlet  or  geometry  effects  (Plee  and  Mellor,  1977) . A recent 
literature  review  of  this  subject  by  Plee  and  Mellor  (1977)  has 
identified  four  possible  mechanisms  which  explain  the  observed  upstream 
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Figure  6-4.  Representation  of  flashback  and  blowoff  limits 

as  a function  of  i)  and  9 for  lean  prevaporized/ 
premixed  combustion  systems. 
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flame  propagation.  All  are  functions  of  combustor  inlet  ten5)erature , 
pressure,  velocity,  local  equivalence  ratio  and  overall  combustor  geom- 
etry; however,  the  latter  two  are  also  related  to  flow  disturbances 
in  the  fuel  preparation  tube  and  even  the  method  of  fuel  injection. 

The  first  mechanism  is  autoignition  which  results  from  the 
gas  residence  time  exceeding  the  fuel  ignition  delay  time.  This  causes 
the  spontaneous  ignition  of  the  fuel/air  mixture  in  the  fuel  preparation 
tube.  Fig.  6-5  illustrates  schematically  three  other  mechanisms  of 
iqjstream  flame  propagation  for  a combustor  utilizing  a sudden  expansion 
flameholder.  Classical  flashback  (as  discussed  in  the  previous  section) 
is  represented  by  Fig.  6- 5a.  In  this  case,  i?)stream  flame  propagation 
usually  occurs  through  the  boundary  layer  and  results  in  a stable  flame 
at  the  fuel  inlet  (Wampler  et  al.,  1976).  Thus  in  classical  flashback, 
the  flame  actually  propagates  upstream,  while  for  autoignition,  sponta- 
neous combustion  occurs.  Unfortunately,  it  is  usually  difficult 
to  distinguish  between  these  two  phenomena. 

Figure  6-5b  and  6- 5c  show  examples  of  flow  disturbances  which 
induce  upstream  combustion.  Flame  stabilization  in  the  separated  flow 
region  can  be  charaterized  by  the  relationship  between  the  characteristic 
dimension  (d)  of  the  flow  disturbance  and  the  distance  (L)  between 
the  disturbance  and  the  flameholder.  When  these  are  the  same  order  of 
magnitude  (d  'v  L) , the  separated  flow  region  interacts  directly  with  the 
flameholder  recirculation  zone  enabling  the  flame  to  propagate  tpstream. 
If  d « L,  the  recirculation  zone  established  by  the  disturbance  does 
not  interact  directly  with  the  flameholder  recirculation  zone.  Instead, 
the  qpstream  zone  provides  sufficient  residence  time  for  the  fuel/air 
mixture  in  the  wake  region  to  autoignite  as  shown  in  Fig.  6- 5c.  This 
differs  from  conventional  autoignition  which  is  a result  of  the  gas 
residence  time  exceeding  the  fuel  ignition  delay  time  in  the  mean  flow. 
Plee  and  Mellor  (1977)  concluded  that  mechanisms  b and  c of  Fig.  6-5 
and  autoignition  were  responsible  for  the  upstream  combustion  reported 
in  the  literature  for  non-catalytic  combustors.  Only  catalytic 
combustors,  which-  operate  at  veiy  low  reference  velocities,  appear 
capable  of  experiencing  classical  flashback. 
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Upstream  flame  propagation  as  a result  of 
flashback  in  the  boundary  layer  or  along 
the  centerline. 
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B.  Upstream  flame  stabilization  as  a result 
of  flame  propagation  through  the  separated 
flow  region  (d  L) . 
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C.  Upstream  flame  stabilization  as  a result  of 
autoignition  in  the  separated  flow  region 
(d  « L) . 


Figure  6-5.  Flame  schematic  illustrating  three  possible 
mechanisms  of  upstream  flame  propagation. 
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Examples  of  other  possible  flow  disturbances  which  cause 
ipstream  flame  propagation  are  shown  schematically  in  Fig.  6-6.  Roffe 
and  Ferri  (1976),  Roffe  (1976)  and  Craig  (1976)  were  able  to  avoid  prob- 
lems of  upstream  flame  stabilization  by  using  a straight  cylindrical 
fuel  preparation  tube  which  eliminated  flow  disturbances  altogether. 
Roberts,?,  et  al.  (1976)  also  experienced  i^stream  combustion  caused  by 
separated  flow  at  some  point  in  the  radial  to  axial  flow  transition  of 
a Vortex  Air  Blast  combustor.  In  fact,  flow  separation  and  finally 
combustion  in  the  fuel  preparation  tube  can  even  be  triggered  by  small 
details  in  the  hardware  such  as  the  failure  to  properly  connect  a slip 
joint  in  the  mixing  tube  (Roffe,  1976).  Only  the  upstream  flame 
propagation  reported  by  Anderson  et  al.  (1975),  Wan^sler  et  al.  (1976), 
Blazowski  and  Bresowar  (1974)  and  Blazowski  and  Walsh  (1975)  in  catalytic 
combustors  are  in  the  region  where  one  would  expect  classical  flashback; 
however,  only  Wampler  et  al.  (1976)  has  completely  eliminated  all  flow 
disturbances . 

In  summary.  Flee  and  Mellor  (1977)  concluded  that  classical 
flashback  has  not  yet  been  observed  in  non-catalytic  combustors  burning 
liquid  jet  fuels  or  propane.  Instead,  upstream  combustion  as  a result 
of  autoignition  and  mechanisms  b and  c of  Fig.  6-5  are  responsible  for 
the  reported  upstream  flame  propagation.  Only  catalytic  combustors  are 
at  a small  enough  0 ^ (see  Fig.  6-4)  for  classical  flashback;  however, 
flashback  and  autoignition  are  usually  indistinguishable.  It  therefore 
appears  that  in  non-catalytic  conibustors  a straight  cylindrical  fuel 
preparation  tube  will  avoid  problems  of  ipstream  combustion  by  elimi- 
nating regions  of  separated  flow.  However,  care  is  advised  since  the 
exact  location  of  the  lean  flashback  limit  has  not  been  determined. 

There  may  be  operating  conditions  at  high  inlet  temperatures  or  pressures 
vdiich  are  within  the  classical  flashback  regime  for  non-catalytic 
combustors.  Further  details  are  given  in  the  review  by  Plee  and  Mellor 
(1977). 

b.  Blowoff 

Unlike  flashback,  a large  quantity  of  experimental  blowoff 
data  currently  exists  for  turbulent  bluff  body  stabilized  flames. 


illustrating  flow  disturbances  which  cause  upstream  flame  stabilization. 
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However,  early  studies  considered  only  sijnple  hydrocarbon  fuels  such  as 
CH^  and  C^Hg  at  relatively  low  inlet  temperatures  and  pressures  compared 
to  practical  combustors.  These  investigations  have  shown  that  in  general 
the  blowoff  limit  widens  with  increasing  pressure  and  ten^ierature  or 
decreasing  reference  velocity;  maximum  flame  stabilization  occurs  near 
stoichiometric  with  approach  flow  equivalence  ratio,  inlet  temperature 
and  reference  velocity  having  the  strongest  effect  (see  discussions  by 
Williams  and  Shipman,  1953;  Longwell,  1953).  Longwell  et  al.  (1949) 
report  a very  weak  influence  of  pressure  on  flame  stabilization.  The 
blowoff  limit  is  also  extended  by  heating  the  flameholder  as  well  as 
decreasing  the  incoming  turbulence  level  (Williams  et  al.,  1949).  In 
addition,  the  inportance  of  approach  flow  turbulence  on  the  limit 
decreases  with  increasing  flameholder  dimension  (Williams  et  al.,  1949; 
Longwell,  1953). 

Flameholder  size,  type  and  shape  also  influence  stabilization. 
In  general,  a two-dimensional  stabilizer  (cylinder  across  a rectangular 
duct)  lias  a much  wider  blowoff  limit  than  a three-dimensional  flameholder 
(circular  disc  in  a duct)  (Williams  and  Shipman,  1953) . Williams  and 
Shipman  (1953)  also  report  that  an  increase  in  the  characteristic 
dimension  of  the  flameholder  leads  to  an  increase  in  the  range  of  stable 
operation;  however,  as  the  size  of  the  stabilizer  becomes  con^arable  to 
the  housing  or  duct  dimension,  a further  increase  narrows  the  blowoff 
limit.  Thus,  an  optimum  flameholder  dimension  exists  for  a given 
combustor  housing.  Flameholder  shape  appears  to  have  only  a second 
order  effect  on  the  extinction  limit  (Longwell  et  al.,  1949;  Williams 
et  al . , 1949 ; Longwell , 1953) . 

Early  correlation  techniques  considered  the  blowoff  velocity 
(Vj^)  as  a function  of  equivalence  ratio  (<(>{jq)  • Variations  in  pressure 
(P),  inlet  ten^ierature  (T^)  and  flameholder  dimension  (D)  were  also 
included  using  equations  of  the  form 
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by  Longwell  et  al.  (1953)  and  DeZubay  (1950).  However,  like  the  velocity 
gradient  versus  equivalence  ratio  correlations  of  the  previous  section  a 
universal  curve  was  not  possible  because  values  of  a,  b,  c,  varied  between 
experimenters . 

Spalding  and  Tall  (1954)  proposed  a Peclet  number  correlation 
similar  to  that  of  Putnam  and  Jensen  (1949)  which  is  derived  by  non- 
dimens ionali  zing  the  energy  equation.  An  example  of  this  correlation 
is  shown  in  Fig.  6-7  where  here  the  Peclet  number  is  based  on  the  free 
stream  velocity  past  the  flameholder  at  blowoff  (Vj^^) , a characteristic 
dimension  (D)  and  the  thermal  diffusivity  (a)  and  is  plotted  against 
the  Peclet  number  based  on  the  laminar  flame  speed.  The  left  most 
curve  in  Fig.  6-7  represents  two-dimensional  flameholders  while  the 
right  most  curve  refers  to  three-dimensional  stabilizers.  The  upper 
part  of  the  axisyiunetric  curve  has  a slope  of  approximately  two  and 
relates  to  the  case  where  turbulent  transport  properties  dominate  the 
stabilization  process;  the  lower  part  of  this  curve  has  a slope  of 
1.5  and  corresponds  to  a dominance  of  molecular  transport.  Variations 
in  pressure,  inlet  temperature  and  fuel  type  are  accounted  for  inplicitly 
through  the  laminar  flame  speed  and  thermal  diffusivity.  However,  as 
discussed  in  the  last  section,  calculation  of  laminar  flame  speeds 
required  for  Pesj^  is  difficult  for  practical  hydrocarbon  fuels  at  elevated 
inlet  temperatures  and  pressures.  Altenkirch  and  Mellor  (1973)  disciiss 
the  difficulty  in  applying  the  Peclet  number  correlation  to  data  obtained 
using  a sinplified  prevaporizing/premixing  combustor  burning  liquid 
propane  at  950K  inlet  temperature  and  4 atm  pressure . 

Spalding  (1955)  has  summarized  various  theories  for  blowoff 
of  bluff  body  stabilized  flames.  In  the  framework  of  characteristic 
times,  two  different  theories  have  evolved.  Longwell  et  al.  (1953) 
state  that  bluff  body  stabilization  occurs  only  if  the  mass  transport 
to  the  recirculation  zone  is  not  great  enough  to  blow  out  an  equiva- 
lent well  stirred  reactor.  Thus,  a stable  flame  exists  if  the  time 
a fluid  particle  spends  in  the  recirculation  zone  is  larger  than  the 
fuel  ignition  delay  time.  The  second  theory  (Zukoski  and  Marble,  1956) 
considers  the  adjacent  shear  layer  rather  than  the  recirculation  zone  as 
tile  region  which  determines  flame  stabilization.  Zukoski  and  Marble  (1956) 
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found  that  the  ten^jerature  in  the  recirculation  zone  is  uniform 
(approximately  nine -tenths  the  adiabatic  flame  teirperature)  and  that 
the  ratio  of  the  wake  tenperature  to  the  adiabatic  flame  teirperature 
remains  approximately  constant  as  the  air  velocity  increases  toward  the 
blowoff  limit.  Furthermore,  the  size  of  the  recirculation  zone  is 
established  by  the  fluid  mechanics  and  is  independent  of  the  detailed 
chemistry;  the  wake  temperature  is  solely  a function  of  the  approacli 
flow  equivalence  ratio  and  inlet  teirperature.  Blowoff  occurs  if  the 
time  a fluid  particle  spends  in  the  shear  layer  is  less  than  the  fuel 
ignition  delay  time. 

Altenkirch  and  Mellor  (1973,  1975)  have  compared  these  two 
theories  in  detail  and  concluded  that  the  shear  layer  controls  blowoff 
in  agreement  with  Zukoski  and  Marble  (1956) . In  addition,  the  Zukoski 
and  Marble  theory  can  explain  the  blowoff- relight  hysterisis  observed  in 
gas  turbine  combustors  (Hottel  et  al.,  1958).  Fig.  6-8  summarizes  the 
physical  model  for  flame  stabilization  proposed  by  Altenkirch  and  Mellor 
(1973,  1975).  Using  the  technique  of  Gosman  et  al.  (1969),  they 
concluded  that  the  time  a fluid  particle  spends  in  the  recirculation 
zone  is  too  long  to  be  characteristic  of  blowoff;  rather,  it  is  the 
purpose  of  the  recirculation  zone  to  provide  heat  and  free  radicals  to 
the  shear  layer  between  the  hot  recirculating  products  and  the  relatively 
cold  unbumed  gases.  If  the  residence  time  in  the  shear  layer  is  less 
than  the  fuel  ignition  delay  time,  the  fresh  mixture  will  never  ignite 
and  the  flame  cannot  stabilize  and  propagate  into  the  free  stream. 

Thus  in  agreement  with  Zukoski  and  Marble  (1956)  the  Damkbhler  simi- 
larity groip  (residence  time/ignition  delay  time)  in  the  shear  layer 
controls  blowoff.  Furthermore  because  the  residence  time  in  the  wake 
region  exceeds  the  shear  layer  residence  time,  the  recirculation  zone 
should  be  extinguished  after  the  shear  layer.  Indeed,  this  has  been 
reported  in  the  early  literature.  Longwell  et  al.  (1949)  and  Williams 
et  al.  (1949)  observed  that  under  certain  operating  conditions  the 
free  stream  flame  blew  off  prior  to  the  wake  region.  Since  the  recir- 
culation zone  represents  only  minor  heat  release  relative  to  the  main 
stream  flame,  this  final  blowout  is  only  of  secondary  interest 
(Blazowski  and  Henderson,  1977). 


Well  Stirred  Reactor,  Long  Residence  Time 
High  Tenperature,  Source  of  Heat  and 
Active  Species 


Physical  Nbdel  of  Flame  Stabilization 
from  Altenkirch  and  Mellor  (1973) . 
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Finally,  Altenkirch  and  Mellor  (1975)  following  Spalding 
(1971)  suggest  that  since  most  flows  of  practical  interest  are  turbulent, 
the  shear  layer  residence  time  should  be  replaced  by  a turbulent  mixing 
time  in  the  shear  layer  region  between  the  recirculating  burned  gases 
and  the  free  stream.  The  characteristic  time  model  therefore  becomes 
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where  t . = turbulent  mixing  time 
and  = ignition  delay  time  of  the  hydrocarbon 

fuel. 

Eq.  6-21  suggests  that  the  hot  turbulent  eddies  present  in  the  shear 
layer  region  must  ignite  before  they  are  quenched  by  the  relatively 
cold  free  stream;  if  they  cannot  ignite,  the  shear  layer  flame  blows 
off.  The  rate  at  which  these  eddies  ignite  is  a function  of  the 

local  fuel/air  ratio  and  ten^erature  and  is  only  related  to  the  chemistry 
of  the  system,  while  the  quenching  rate  is  controlled  solely  by 

the  fluid  mechanics. 

Following  Tuttle  et  al.  (1975,  1976,  1977),  the  turbulent 

mixing  time  is  assumed  to  be  the  characteristic  breakdown  time  of  the 

large  scale  eddies,  is  a function  of  an  integral  length  scale 

divided  by  the  turbulent  (rms)  velocity  evaluated  in  the  appropriate 

region  of  interest  (see  Chapter  II).  The  integral  length  scale  and 

turbulent  velocity  generally  scale  with  a geometric  dimension  in  the 

combustor  and  a convective  velocity,  respectively.  Unfortunately, 

appropriate  geometric  length  scales  cannot  be  specified  a priori 

but  must  instead  be  determined  through  analysis  of  the  experimental 

data  (Tuttle  et  al.,  1976).  For  blowoff,  the  turbulent  mixing  time 

is  taken  to  be  a ratio  of  the  flameholder  width  (L) , which  is  also  a 

measure  of  the  size  of  the  recirculation  zone,  divided  by  the  air 

velocity  at  the  edge  of  the  stabilizer  (V  ™) , 

a,  1 


V 


a,T 


(msec) 


6-22 
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(prime  denotes  tenperature  dependence  in  velocity) , since  according  to 
Zukoski  and  Marble  (1956)  this  is  the  region  where  stabilization  occurs. 
The  annulus  velocity  is  confuted  from  the  continuity  equation. 
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and  is  a function  of  the  reference  velocity  blockage  ratio 

(B) . In  addition,  the  velocity  is  evaluated  at  the  tenqjerature  in 
the  shear  layer  (T)  which  must  be  higher  than  the  inlet  temperature 
because  it  is  adjacent  to  a hot  recirculation  zone.  Zukoski  and  Marble 
(1956)  have  shown  that  the  recirculation  zone  maintains  its  tenqierature 
regardless  of  the  velocity  past  the  flameholder  even  up  to  blowoff. 

Eq.  6-22  represents  the  same  shear  layer  mixing  time  that  Tuttle  et  al. 
(1977)  used  for  the  correlation  of  carbon  monoxide  emissions  except 
here  the  inner  part  of  the  shear  layer  rather  than  the  cold  outer  region 
is  considered.  Therefore,  the  ratio  T/T.  * in  Eq.  6-23  accounts  for 
the  acceleration  of  the  air  flow  in  the  shear  layer  region  as  a result 
of  the  increased  temperature. 

The  ignition  delay  time  for  a lean  reacting  eddy  is 


.E/RT 


^hc 


(sec) 
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from  the  previous  laminar  flow  analysis.  Combining  Eq.  6-22  and  6-16, 
the  characteristic  time  model  for  lean  blowoff  becomes 


6-24 


Now  since  V is  a function  of  the  recirculation  zone  tenperature 
(approach  flow  equivalence  ratio) , we  write  the  ratio  T/T^  on  the 
right  hand  side  of  Eq.  6-24  so  that  the  fluid  mechanics  and  chemistry 

Tf-y 

Since  velocities  here  are  much  higher  than  for  laminar  flow,  increased 
coii?)ressibility  effects  must  also  be  included. 
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are  again  uncon)led, 
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For  the  prevaporizing/premixing  combustor  shown  in  Fig.  6-8, 
the  appropriate  length  in  should  scale  with  the  housing  diameter 
(op  minus  the  baffle  diameter  (D)  since  this  is  a measure  of  the  size 
of  the  recirculation  zone.  The  temperature  in  is  assumed  propor- 
tional to  the  adiabatic  flame  temperature  at  the  approach  flow 
equivalence  ratio  just  prior  to  blowoff  (an  estimate  of  the  thermal 
driving  force  between  the  recirculation  zone  and  the  shear  layer) 
following  Zukoski  and  Marble  (1956) . The  model  can  also  be  extended 
to  three  other  simplified  flameholder  configurations  shown  schematically 
in  Fig.  6-9.  Each  is  characterized  by  a shear  layer  flame  which  is 
ignited  through  transfer  of  heat  and  mass  from  a central  recirculation 
zone  (boundary  shown  with  dashed  lines) . 

Ballal  and  Lefebvre  (1977)  have  recently  obtained  lean 
blowoff  data  from  bluff  body  stabilized  flames  operating  at  inlet 
ten^eratures  and  velocities  typical  of  continuous  combustors  (Fig.  6-9A) ; 
these  conditions  are  summarized  in  Table  6-2.  In  this  study  gaseous 
propane,  introduced  upstream  of  the  conical  baffle,  is  mixed  with  air 
prior  to  combustion  eliminating  any  heterogeneous  effects.  Variations 
in  blowoff  equivalence  ratio  with  pressure,  inlet  ten^erature,  reference 
velocity,  baffle  size,  housing  diameter  and  blockage  ratio  have  been 
examined  independently.  The  lean  limit  was  obtained  by  holding  all 
inlet  conditions  constant  and  slowly  decreasing  the  fuel  flow  until 
blowoff  was  observed  visually. 


Configuration  A Ballal  and  Lefebvre  (1977) 


^^__recirculat 

— 

~'S  boundary 

air 

V_  shear  layer 

1 XCUIl^ 


Configuration  B AFOSR  Combustor  (Fig.  3-3) 


Configuration  C Plee  and  Mellor  (1978a) 

(negligible  heterogeneous  effects) 


Figure  6-9.  Schematic  of  three  simplified  bluff 
body  stabilized  flames. 
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Table  6-2.  Burner  Operating  Conditions 
(Ballal  and  Lefebvre,  1977) 


Reference  Velocity 
Pressure 

Inlet  Temperature 
Pipe  Diameter 
Flameholder  Blockage  Ratio 
Fuel 

Geometry 

Incoming  Turbulence  Intensity 


10  - 100  m/sec 
0.2  - 0.9  atm 
300  - 575K 
10,  15,  25  cm 
4 - 34% 

CjHg  (gas) 

45°  conical  baffle 
up  to  15% 


Configuration  B represents  the  present  AFOSR  geometry 
discussed  in  Chapter  III  and  configuration  C is  the  disc-in-duct 
flameholder  studied  extensively  by  Tuttle  et  al.  (1975,  1976,  1977) 
for  EPA.  Plee  and  Mellor  (1978a)  have  recently  obtained  lean  blowoff 
data  for  configuration  C which  consists  of  a disc  mounted  in  a duct  with 
liquid  fuel  injected  at  the  center  via  a simplex  pressure  atomizing 
nozzle.  The  disc- in-duct  combustor  is  fundamentally  different  from 
these  discussed  previously  because  the  fuel  is  not  prevaporized  and 
premixed  upstream  of  the  flameholder. 

Appropriate  characteristic  time  scaling  parameters  for  each 
of  these  combustors  are  listed  in  Table  6-3.  The  turbulent  mixing  time 
which  is  a function  of  L and  is  relatively  siti^Jle  to  determine  for 
each  of  these  flameholder  configurations.  In  each  case,  is  the  ratio 
of  the  flameholder  width  divided  by  the  annular  velocity  (V  ) . Thus, 
the  appropriate  length  scale  for  the  tube-and-disc  con±)ustor  is  assumed 
to  be  the  baffle  diameter  (D)  minus  the  tube  diameter  (d) . Note  that  the 
quantity  (D-d)  which  is  related  to  the  size  of  the  recirculation  zone 
goes  to  the  proper  limit  (D)  as  the  tube  diameter  approaches  zero. 

Parameters  relevant  to  the  ignition  delay  time  are  not  always 
as  straightforward  to  evaluate.  Since  the  combustor  of  Fig.  6-9A  is 
essentially  identical  to  the  flameholder  shown  in  Fig.  6-8  except  for  the 
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location  of  the  bluff  body,  the  appropriate  equivalence  ratio  and  flame 
temperature  are  evaluated  in  the  same  manner.  However,  for  configura- 
tion B,  we  postulate  that  the  inner  shear  layer  (shear  layer  beginning 
at  the  edge  of  tlie  tube.  Fig,  6-9B)  is  responsible  for  flame  stabiliza- 
tion since  high  fuel  concentrations  are  present  at  this  point.  Hot 
recirculating  burned  gases  should  also  maintain  this  region  at  a 
relatively  high  temperature  based  on  the  ten^jerature  measurements  of 
Chapter  IV.  Therefore,  application  of  Eq.  6-25  requires  an  appropriate 
length  scale,  velocity,  equivalence  ratio  and  ten^erature  evaluated  in 
this  regime.  The  equivalence  ratio  representative  of  the  tube-and-disc 
combustor  is  assumed  to  be  (})^  which  is  based  on  the  air  flow  split 
through  the  tube  rather  than  the  total  air  flow  rate  (4*^^)  • is 
richer  than  the  overall  equivalence  ratio  and  is  a better  measure  of 
the  fuel  concentration  near  the  shear  layer  adjacent  to  the  tube  exit. 

By  analogy  with  the  blowoff  model  of  the  previous  section, 
the  appropriate  temperature  for  should  be  the  adiabatic  flame 
ten^ierature  (T(j,^)  based  on  the  equivalence  ratio  in  the  fuel  prepara- 
tion tube.  However,  this  does  not  take  into  account  the  influence  of 
the  relatively  cold  inlet  air  flowing  around  the  outside  of  the  disc. 
Obviously  there  exists  a fundamental  difference  between  the  configura- 
tions studied  by  Ballal  and  Lefebvre  (1977)  and  the  tube-and-disc 
flameholders.  In  the  latter  case,  the  recirculation  zone  is  not  com- 
pletely surrounded  by  a hot  reacting  shear  layer;  instead,  the  outer 
shear  layer  at  the  edge  of  the  disc  is  almost  entirely  free  of  fuel  and 
should  be  relatively  unreactive.  Therefore,  a better  representation  of 
the  tenperature  of  the  inner  shear  layer  is  an  average  between  T,j,^  and 
the  inlet  temperature.  This  simplified  analysis  allows  for  the  quenching 
effect  of  the  cold  outer  recirculation  zone  (see  Fig.  6-9B).  If  the 
disc  diameter  equals  the  housing  diameter,  then  the  entire  recirculation 
zone  would  be  surrounded  by  either  a wall  or  a reacting  shear  layer  (the 


dunp  ramjet  combustor)  and  the  appropriate  ten^erature  would  then  be  T(j,^ 
(equal  to  T,j,^y)  as  in  the  case  of  Fig.  6-8. 

Configuration  C differs  substantially  from  the  previous 
flameholders  studied  because  the  fuel  is  injected  at  the  center  of  the 
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recirculation  zone  and  is  not  therefore  prevaporized  or  premixed  upstream. 
Tuttle  et  al.  (1976,  1977)  successfully  treated  this  configuration  as 
a turbulent  diffusion  flame  for  modeling  NO^  formation.  Their  results 
show  that  very  high  fuel/air  ratios  exist  in  the  recirculation  zone  due 
to  fuel  vaporization  and  then  decrease  into  the  main  stream.  Maximum 
reaction  rates  and  very  high  tenperatures  characteristic  of  NO  formation 
are  created  by  the  stoichiometric  contour  present  in  the  shear  layer 
region  between  the  fuel  rich  rec  .^rculation  zone  and  the  free  stream. 

These  stoichiometric  regions  also  account  for  the  observations  that 
conventional  diffusion  flame  combustors  are  generally  more  stable  than 
prevaporizing  premixing  combustors.  Therefore  since  the  maximum 
reaction  rate  occurs  primarily  around  the  stoichiometric  contour,  the 
appropriate  flame  temperature  in  should  be  evaluated  at  <()=1. 

The  equivalence  ratio  in  is  not  evaluated  at  ((>=1.0  however  since 
all  of  the  eddies  are  not  at  stoichiometric.  Instead,  the  percent 
of  stoichiometric  eddies  controlling  flame  stabilization  should  be 
proportional  to  the  overall  equivalence  ratio.  Thus  as  (p  increases, 
the  ignition  delay  time  becomes  smaller  and  the  flame  becomes  more 
stable.  Examples  of  this  particular  correlation  are  shown  in  Plee  and 
Mellor  (1978b). 

Application  of  the  characteristic  time  model  and  evaluation 
of  for  blowoff  are  demonstrated  first  with  the  flameholder  of 
Fig.  6-9A.  Using  the  experimental  data  of  Ballal  and  Lefebvre  (1977), 
the  ignition  delay  time  becomes 

, „ 21000/RT 

Tj*  = 10  ^ s; — T (msec)  6-26 

A 

where  the  prime  indicates  that  contains  the  ratio  T/Tj^  from  the 
annulus  velocity.  Note  that  both  the  pre-exponential  (10"^)  and  the 
activation  energy  (E  = 21,000  cal/mole)  are  determined  enpirically. 

An  activation  energy  of  21,000  cal/mole  gives  the  best  correlation 

It 

Failure  to  include  the  ratio  T/Tin  from  the  velocity  term  decreases 
the  correlation  coefficient  in  Fig.  6-10  from  0.98  to  0.93. 
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and  compares  favorably  with  values  of  24,600  cal/mole  and  24,200  cal/mole 
reported  by  Engleman  et  al.  (1973)  and  Edelman  (1977),  respectively. 

Figure  6-10  illustrates  the  characteristic  time  model  for  the 
lean  blowoff  data  of  Ballal  and  Lefebvre  (1977) . Here  the  shear  layer 
mixing  time  and  the  fuel  ignition  delay  time  are  evaluated  using  Eq.  6-25 
with  evaluated  from  Eq.  6-26.  The  temperature  in  Eq.  6-25  and  6-26 
is  taken  to  be  the  adiabatic  flame  temperature  at  the  approach  flow 
equivalence  ratio  just  prior  to  blowoff  (Table  6-3)  and  is  conputed 
using  a version  of  the  equilibrium  computer  program  of  Gordon  and 
McBride  (1971).  The  correlation  is  shown  to  be  linear  with  a y-intercept 
approximately  zero  as  the  theory  predicts.  Thus,  one  approaches  the 
lean  blowoff  limit  by  either  increasing  the  annular  velocity  or 
decreasing  the  flameholder  characteristic  dimension  (fT^jj^) , or 
decreasing  the  overall  equivalence  ratio  • 

Analysis  of  the  turbulent  mixing  time  predicts  an  optimum 
disc  diameter  for  maximum  stability  at  a given  power  setting  which 
is  solely  a function  of  the  housing  diameter.  Since 


r = 2_ 
si  V 


D(D^  - D^) 
c 

7~ 
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maximum  flame  stabilization  exists  is  a maximum)  at 
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Thus,  increasing  the  characteristic  dimension  of  the  flameholder  leads 

to  an  increase  in  the  range  of  stable  operation;  however,  as  the  size 

of  the  stabilizer  exceeds  B//T,  a further  increase  narrows  the  blowoff 

c 

limit  (■*’Tgj^)  as  reported  by  Williams  and  Shipman  (1953).  Eq.  6-28 
is  also  identical  to  the  equation  obtained  by  Tuttle  et  al.  (1977) 
for  minimum  CX)  emissions. 

In  sunmary,  a characteristic  time  model  has  been  developed 
based  on  the  flame  stabilization  theory  of  Zukoski  and  Marble  (1956). 

for  lean  blowoff  is  identical  to  the  mixing  time  obtained  by 
Tuttle  et  al.  (1977)  for  CD  emissions;  however,  kinetic  times  are  quite 


data  by  Ballal  and  Lefebvre  (1977') 


Figure  6-10.  Characteristic  time  correlation  for  lean 
blowoff  in  a bluff  body  stabilized  flame 
burning  gaseous  propane. 
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different.  We  now  seek  to  apply  this  model  to  the  blowoff  data  obtained 
on  the  tube-and-disc  combustor  after  a brief  discussion  of  flashback. 

B.  Experimental  Results  and  Discussion 

1.  Flashback 

For  the  simplified  prevaporizing/premixing  combustors  of 
Fig.  3-3  and  Table  3-2,  neither  flashback  nor  autoignition  were 
observed  using  the  burner  operating  conditions  of  Table  3-4.  Spadaccini 
(1976)  has  shown  that  the  ignition  delay  time  of  JP-4  under  these 
conditions  is  greater  than  10  msec  and  therefore  one  would  not  expect 
autoignition  to  be  a problem  (see  Table  3-2).  Furthermore,  the  absence 
of  classical  flashback  is  consistent  with  the  recent  review  by  Flee 
and  Mellor  (1977)  summarized  in  Section  A of  this  Chapter.  Classical 
flashback  should  not  be  a problem  in  non-catalytic  combustors  utilizing 
straight  cylindrical  fuel  preparation  tubes.  This  means  that  for  the 
combustors  examined  in  this  facility  the  minimum  6 ^ has  not  yet  been 
reached  (see  Fig.  6-4);  therefore,  increasing  the  equivalence  ratio  in 
the  tube  through  stoichiometric  will  not  induce  classical  flashback. 

Although  the  flashback  limit  has  not  been  obtained  to  date, 
these  results  are  important  because  they  help  define  regions  above  this 
limit.  The  apparent  overlap  between  these  stable  operating  conditions 
and  the  reported  "flashbacks"  in  the  literature  is  attributed  to  flow 
disturbances  and  autoignition  rather  than  classical  flashback  (Flee 
and  Mellor,  1977).  Thus,  elimination  of  flow  disturbances  in  the 
fuel  preparation  tube  avoids  problems  of  upstream  combustion  altogether. 

2.  Blowoff 

Blowoff  data  obtained  for  the  five  prevaporizing/premixing 
combustors  of  Fig.  3-3  and  Table  3-2  under  the  operating  conditions 
listed  in  Table  3-4  are  presented  in  Fig.  6-11  through  6-15.  Here 
the  equivalence  ratio  at  blowoff  is  plotted  versus  the  inverse 
theta  parameter  defined  by  Eq.  6-19.  In  addition,  effects  of  vitia- 
tion are  included  using  an  equation  defined  by  Altenkirch  and  Mellor 
(1975) 
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where  dv  = blowoff  equivalence  ratio  including 

vitiation 

<)>^^  = apparent  blowoff  equivalence  ratio 

and  (}»u  ® overall  equivalence  ratio  in  slave 

burner . 

Eq.  6 -29  corrects  the  blowoff  equivalence  ratio  for  vitiation  assuming 
complete  combustion  in  the  preheat  burner. 

Blowoff  equivalence  ratios  obtained  burning  liquid  propane  in 
geometry  A (base  geometry)  are  presented  in  Fig.  6-11.  This  plot 
illustrates  that  the  inverse  theta  pargimeter  correlates  the  data  reason 
ably  well  for  varitaions  in  pressure,  inlet  tenperature  and  velocity. 

A similar  plot  is  available  for  the  same  flameholder  configuration 
burning  Jet  A (Fig.  6-12)  except  here  variations  in  initial  fuel  spray 
S^D  are  also  included  and  are  not  collapsed  by  9 An  increase  in  the 
size  of  the  injector  leads  to  a reduction  in  the  stable  operating  range 
because  the  initial  fuel  spray  SMD  also  increases.  The  propane  and 
Jet  A blowoff  data  are  plotted  together  in  Fig.  6-13  at  one  particular 
inlet  temperature  to  further  show  the  importance  of  atomization  on  the 
limit.  For  liquid  propane,  the  fuel  drops  are  expected  to  flash 
vaporize  iqjon  injection  under  the  burner  inlet  conditions  of  Table  3-4 
(SMD  'V'  0 at  tube  exit) . Jet  A on  the  other  hand  is  not  expected  to 
conpletely  vaporize  prior  to  combustion  as  discussed  in  Chapter  IV,  and 
therefore,  increased  fuel  volatility  also  influences  the  degree  of 
pre vaporization . 

According  to  Fig.  6-13,  heterogeneous  effects  are  only 
important  relative  to  blowoff  at  large  values  of  0 ^ (high  velocity) . 
This  is  attributed  to  the  fact  that  at  high  air  velocities,  less 
residence  time  is  available  in  the  mixing  tube  for  vaporization, 
while  at  low  velocities,  the  residence  time  is  sufficient  such  that 
vaporization  does  not  influence  the  blowoff  limit.  Therefore,  the 
more  fuel  vaporized  in  the  fuel  preparation  tube  prior  to  ignition, 
the  higher  the  local  fuel  to  air  ratio  and  the  wider  the  blowoff  limit. 
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Figure  6-11.  Blowoff  equivalence  ratio  versus  9 
(Propane,  Geometry  A) . 


Figure  6-12.  Blowoff  equivalence  ratio  versus  9 
(Jet  A,  Geometry  A) . 


Figure  6-13.  Blowoff  equivalence  ratio  versus  9 
fPropane  and  Jet  A,  Geometry  A). 
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In  summary,  both  changes  in  fuel  type  and  injector  size  can  have  an 
adverse  effect  on  the  blowoff  limit  if  the  initial  fuel  spray  SMD 
increases . 

Geometry  v^ariations  are  shown  in  Fig.  6-14  and  6-15  for 
combustors  burning  propane  and  Jet  A,  respectively.  In  both  figures, 
injector  size  and  inlet  temperature  are  held  constant  in  order  to  better 
illustrate  geometry  changes . Although  the  inverse  theta  parameter  also 
does  not  correlate  geometry  effects,  it  does  illustrate  certain  trends. 
First  of  all,  the  small  tube  flameholders  (geometries  D and  E)  are  more 
stable  tlian  the  4.4  cm  diameter  tube  configurations  (geometries  A,  B 
and  C)  for  both  fuels.  Also,  decreasing  the  disc  diameter  for  a given 
mixing  tube  diameter  enhances  flame  stabilization  (geometry  B versus 
geometries  A and  C;  geometry  D versus  geometry  E) . Another  interest- 
ing result  is  that  for  combustors  burning  Jet  A the  shorter  tube 
(geometry  C)  has  the  same  apparent  limit  as  the  base  geometry.  This 
can  be  explained  from  examination  of  the  characteristic  droplet 
evaporation  time  introduced  in  Chapter  II. 

For  the  tube-and-disc  combustor,  Eq.  2-11  must  be  modified 

to  account  for  prevaporization  occurring  in  the  fuel  preparation  tube. 

Using  Eq.  2-10,  the  droplet  diameter  at  the  exit  of  the  prevaporizing 

tube  (d  ) becomes 
e 


<“o  ■ Vs 


6-30 


where  B.  = evaporation  coefficient  evaluated  under 
ambient  conditions  in  the  tube 

and  = droplet  residence  time  in  tube. 

Since  droplet  sizes  of  d rather  than  d^  are  exposed  to  the  combustion 

environment,  Eq.  2-11  now  becomes 


d^  - B.  T 
_ o t res 

'^eb  B 
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Thus  if  the  fuel  (propane)  is  conpletely  vaporized  in  the  fuel  prepara- 
tion tube,  the  numerator  of  Eq.  6-31  is  zero  and  can  be  neglected  in 
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14.  Blowoff  equivalence  ratio  versus  9 
(Propane,  all  geometry  variations). 


Blowoff  equivalence  ratio  versus  9 
(Jet  A,  all  geometry  variations^ . 
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tlie  combustion  correlations.  However,  if  conplete  vaporization  does 
not  occur  in  the  mixing  tube,  then  some  fuel  droplets  will  leave  the 
tube  exit  plane  and  enter  the  stabilization  zone  shear  layer.  Further- 
more, the  temperature  in  the  shear  layer  should  be  much  higher  than 

inside  the  tube  so  that  if  both  the  ratio  g^/6  and  t are  small, 

t c res  ’ 
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as  if  the  prevaporization  tube  had  no  influence  whatsoever.  Eq.  6-32 
applies  only  to  circumstances  when  extremely  poor  vaporization  exists 
in  the  tube  and  explains  why  further  shortening  of  the  tube  can  have 
a negligible  effect  on  the  blowoff  limit. 

Having  presented  the  data  using  conventional  gas  turbine 
correlating  techniques,  we  now  endeavor  to  develop  the  characteristic 
time  model  to  better  explain  quantitatively  the  observed  influence  of 
variations  in  geometry  and  atomization  as  well  as  inlet  conditions. 

The  length  scale  for  the  turbulent  mixing  time  is  assumed  to 

be  the  disc  diameter  (D)  minus  the  tube  diameter  (d)  since  the  flame - 
holder  size  (D)  successfully  correlated  data  from  Ballal  and 
Lefebvre  (1977).  Following  Table  6-3,  the  appropriate  velocity  should 
be  the  annular  velocity  (velocity  at  the  tube  exit)  which  for  this 
particular  configuration  can  easily  be  computed  from  the  flameholder 
blockage  (see  the  pitot  tube  measurements  of  Qiapter  IV) . Therefore , 
the  shear  layer  mixing  time  is 


D - d 


si 
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According  to  Table  6-3  the  appropriate  temperature  and 
equivalence  ratio  are  (T(j,^  + T^)/2  and  (j)^,  respectively.  Thus,  the 
characteristic  time  model  for  the  tube-and-disc  becomes 


D - d „ T 

-y 'X-  10  ^ 

a in 


,21000/RT 
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where  <J>  = 

and  T = (T|^  + T.^)/2  . 

Vitiation  effects  are  omitted  when  calculating  the  tube 
equivalence  ratio  because  the  original  derivation  contains  the 
mole  fraction  of  fuel,  rather  than  (f).  The  value  of  remains  unchanged 
with  vitiation  (same  number  of  moles  of  fuel  are  present;  only  the  conposi- 
tion  of  the  air  changes)  and  is  related  to  either  of  the  following 
proportionalities : 


Xjj^  'V  fuel/air  'v  (f/a)^  (|) 


) 


apparent 
vitiated 
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Since  the  stoichiometric  fuel/air  ratio  for  (tiapparent  approximately 
constant  with  typical  hydrocarbon  fuels  while  f/a]^  changes  with 
inlet  temperature,  Eq.  6-35  simplifies  to 


“^app  ^ “^actual 


f(T.  ) 
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Both  aie  valid  and  equivalent  but  obviously  X^^^  <l>gpp  is  much  sinpler 
to  use.  Vitiation  effects  must  however,  be  included  in  flame  tenpera- 
ture  calculations  since  the  presence  of  water  vapor  and  carbon  dioxide 


in  the  approach  flow  decreases  the  near  stoichiometric  flame  temperature 
substantially  at  inlet  temperatures  of  BOOK  (Blazowski  and  Henderson, 
1977). 


In  general,  tube  equivalence  ratios  are  fuel  lean;  however,  under  some 
circumstances,  exceeds  unity  at  blowoff.  Thus,  the  stabilization 
curve  will  have  the  same  shape  as  the  laminar  flow  blowoff  curve  in 
air  (see  Fig.  6-1)  due  to  mixing  of  the  O2  rich  outer  flow  with  the 
fuel  rich  mixture  from  the  tube.  This  form  is  further  evidenced  by 
the  ability  to  obtain  combustion  efficiency  data  at  tube  equivalence 
ratios  in  excess  of  2.0  (Table  3-4).  When  blowoff  occurs  at  greater 
than  one,  the  stoichiometric  flame  temperature  replaces  Ta^  in  Eq.  6-34 
to  account  for  the  mixing  of  oxygen  and  subsequent  reduction  in 
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However,  it  has  also  been  observed  (Fig.  6-13)  that 
atomization  influences  flame  stabilization.  Given  the  same  mixing 
time  as  the  purely  homogeneous  sustem,  the  heterogeneous  combustion 
process  must  both  evaporate  the  fuel  and  ignite  the  fuel  vapor  in  the 
same  time  span.  Moses  (1975)  has  also  reported  that  poor  atomization 
has  a detrimental  effect  on  blowoff.  Thus,  droplets  represent  a 
perturbation  on  the  flame  stabilization  process  which  narrow  the  limit. 
To  account  for  heterogeneous  effects,  Eq.  6-24  is  modified  to 


T „ '''Tv 

Si  he 


+ ax 


eb 


6-37 


2 

where  x^j^  is  conputed  from  the  "d  law"  discussed  in  Chapter  II  and 
evaluated  under  conditions  in  the  shear  layer.  The  parameter  "a"  in 
Eq.  6-37  is  chosen  empirically  to  give  the  best  fit  to  the  experimental 
data.  Finally,  moving  the  ratio  T/T^^  to  the  right  hand  side  as  before 
yields 


"si 


■^hc 


6-38 


The  experimental  data  obtained  for  the  base  geometry  burning 
both  CjHg  and  Jet  A are  shown  in  Fig.  6-16  using  the  model  of  Eq.  6-38. 
The  correlation  includes  variations  in  pressure,  inlet  temperature, 
velocity,  fuel  type  and  injector  size.  It  can  be  seen  that  one 
approaches  the  blowoff  limit  by  increasing  the  velocity  (+  x^^^) , 
increasing  the  initial  droplet  size  (f  x^j^)  and  decreasing  the  tube 
equivalence  ratio  (1  Tj!^^).  Tlie  model  also  predicts  that  droplet 
effects  will  become  inportant  at  high  velocities  (small  x^^^)  since  at 
this  point  Xj^^  and  x^j^  are  the  same  order  of  magnitude.  This  is 
verified  experimentally  in  Fig.  6-13;  droplet  effects  are  most 
significant  at  large  values  of  0 ^ which  correspond  to  high  air 
velocities . 

Geometry  variations  are  shown  in  Fig.  6-17  and  6-18  for 
propane  and  Jet  A,  respectively.  Again  the  correlation  coefficient 
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Figure  6-16.  Characteristic  time  correlation  for  blowoff 
(vJet  A and  Propane , Geometry  A) . 
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is  acceptable  with  the  y- intercept  passing  through  zero  to  within  one 
standard  deviation.  All  of  these  data  are  presented  in  Fig.  6-19  and 
correlated  with  the  characteristic  time  equation 


T „ = 2.24  (t'  + O.OIIt'.  ) + 0.06 

s£  he  eb-^ 


6- 


vdiich  represents  the  blowoff  limit.  Knowledge  of  the  flameholder 
configuration,  inlet  conditions,  and  atomization  determines  and 
from  these,  can  easily  be  cojrputed  which  is  a function  of 
the  inlet  conditions  and  the  blowoff  equivalence  ratio.  The  blowoff 
velocity  is  also  calculated  in  a similar  fashion  if  the  geometry, 
pressure,  inlet  temperature,  atomization  and  equivalence  ratio  are 
given.  In  addition,  this  equation  can  predict  simple  flameholder 
modifications  necessary  to  optimize  flame  stabilization  for  a given 
operating  condition.  The  cliaracteristic  time  model,  although  not  as 
straightforward  as  conventional  loading  parameters,  does  include 
variations  in  geometry  and  atomization  as  well  as  inlet  conditions 
with  a minimum  amount  of  algebra. 


3.  Shrouded  Combustor 

The  lean  blowoff  limit  has  also  been  observed  for  the 
shrouded  combustor  of  Fig.  3-4  operating  at  inlet  conditions  described 
in  Table  3-5.  Table  6-4  summarizes  these  results  which  seem  to 
indicate  that  the  shrouded  combustor  is  more  stable  than  the  simplified 
version  (vitiation  effects  have  been  included) . This  apparent  decrease 
in  the  blowoff  equivalence  ratio  is  probably  a result  of  the  stream- 
lining effect  of  the  shroud  as  indicated  by  the  pitot  measurements  of 
Chapter  IV.  However,  it  is  only  of  secondary  inportance  under  these 
burner  operating  conditions.  The  limit  is  still  primarily  influenced 
by  flameholder  size,  reference  velocity  and  inlet  temperature. 
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Figure  6-19.  Characteristic  time  correlation 
for  blowoff  (all  data) . 
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Table  6-4.  Effect  of  Shrouded  Conbustor  on  Blowoff 
(Jet  A,  Geometry  A) 


P(atm) 

T.  (K) 
in  ^ 

Vref(n>/sec) 

'^bo.vit 

Con±)ustor 

8 

800 

17.0 

0.18  - 0.21 
0.15  - 0.16 

sin5)lified 

shrouded 

4 

800 

34.0 

0.28  - 0.30 
0.25  - 0.26 

simplified 

shrouded 

C.  Summary 

The  characteristic  time  model  has  been  extended  to  include 
blovroff , first  for  a conical  flameholder  burning  gaseous  propane 
(Ballal  and  Lefebvre,  1977)  and  later  for  the  AFOSR  tube-and-disc 
configuration.  Important  parameters  in  the  characteristic  time  model 
of  Eq.  6-25  for  both  combustors  are  summarized  in  Table  6-3. 

These  data  are  plotted  together  in  Fig.  6-20  using  the 
appropriate  parameters  of  Table  6-3.  Agreement  between  the  two  sets 
of  data  is  remarkably  good.  The  characteristic  time  model  correlates 
variations  in  pressure,  inlet  temperature,  velocity,  fuel  type, 
injector  and  geometry  for  the  two  different  configurations  using  the 
equation 


= 1.88  (Xj^^  + O.Ollx^j^)  + 0.07  6-40 

(prime  denotes  that*the  ratio  T/T^^  in  the  velocity  term  has  been 
included  on  the  right  hand  side  of  the  equation),  x . in  Eq.  6-40  is 
identical  to  the  mixing  time  derived  by  Tuttle  et  al.  (1977)  for  CD 
emissions;  however,  the  ignition  delay  time  is  quite  different  from 
the  characteristic  time  for  CO  oxidation.  In  addition  to  the  data 
presented  here,  blowoff  limits  recently  obtained  by  Plee  and  Mellor 
(1978a)  for  a disc  stabilized  flame  with  negligible  heterogeneous 
effects  can  also  be  included  in  the  correlation  (Plee  and  Mellor,  1978b). 
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Figure  6-20.  Characteristic  time  model  for  lean 
blovroff  data  obtained  from  bluff 
body  stabilized  flames. 
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Flameholder  shape  is  not  included  in  the  model  since  results  in  the 
previous  section  indicate  that  it  is  essentially  negligible;  our  ability 
to  correlate  these  two  different  configurations  enforces  this  belief. 

We  must  mention,  however,  that  it  is  not  always  a simple 
procedure  to  obtain  the  correct  parameters  for  the  characteristic  time 
model.  The  key  to  success  is  the  proper  identification  of  important 
regions  of  the  flow  field  either  from  examination  of  the  literature  or 
from  detailed  probing.  Exan^les  of  some  variations  in  the  model  which 
do  not  correlate  as  well  are  shown  in  Table  6-5.  Of  all  the  correlations 
examined  (Table  6-5  lists  only  a few).  Fig.  6-19  gives  the  best  correlation 
coefficient  (r)  and  also  has  approximately  the  same  slope  (m)  as  the  data 
of  Ballal  and  Lefebvre  (1977).  Therefore,  the  characteristic  time  model 
is  useful  to  the  designer  for  not  only  correlating  data  but  also  possibly 
identifying  key  regions  of  the  flow  field. 

Since  the  flashback  mechanism  is  very  similar  to  blowoff , the 
characteristic  time  model  sliould  also  be  applicable  except  that  for 
flashback,  times  are  evaluated  inside  the  fuel  preparation  tube  rather 
than  the  recirculation  zone  shear  layer.  However,  examination  of  current 
literature  indicates  that  classical  flashback  has  not  been  observed  in 
non-catalytic  combustors  burning  hydrocarbon  fuels  but  rather  upstream 
flame  propagation  is  a result  of  autoignition  and  flow  disturbances  in 
the  mixing  tube.  In  agreement  with  the  literature,  flashback  has  not 
been  observed  in  this  study  using  a straight  cylindrical  fuel  preparation 
tube. 


Table  6-5.  Conparison  of  Various  Parameters  in  the  Model 
(tube-and-disc  geometry,  E = 21,000  cal/mole) 
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VII.  FUTURE  EFFORTS 

One  primary  objective  of  the  program  continuation  is  the 
relocation  of  the  Purdue  combustion  tunnel  to  an  adjacent  building  in 
order  to  take  advantage  of  an  indirect  fired  air  heater  and  con^uterized 
data  acquisition  system.  These  major  inprovements  should  greatly 
increase  the  quantity  as  well  as  the  reliability  and  accuracy  of  our 
experiments . Assuming  that  the  new  laboratoiy  is  operational  in  the 
near  future,  detailed  experiments  to  better  characterize  the  combus- 
tion flow  field  and  flame  stabilization  process  are  proposed.  Lean 
blowoff  limits  have  already  been  successfully  correlated  using  the 
characteristic  time  model;  however,  the  classical  flashback  limit 
still  remains  unknown  for  our  prevapori zing/premixing  combustors.  This 
is  due  to  the  fact  that  the  minimum  0 ^ (velocity)  attainable  in  the 
current  facility  is  limited  by  flame  stabilization  in  the  pre-heat 
combustor.  The  indirect  fired  air  heater  in  the  new  laboratory  not 
only  eliminates  vitiation  effects  but  also  allows  much  lower  reference 
velocities  needed  for  determination  of  the  classical  flashback  limit. 

In  addition,  results  of  the  combustion  efficiency  section 
clearly  indicate  the  need  for  a better  understanding  of  the  combustion 
flow  field.  To  date,  the  characteristic  time  model  for  combustion 
efficiency  cannot  include  variations  in  pressure,  inlet  temperature, 
equivalence  ratio,  air  velocity,  fuel  type  and  geometry  simultaneously 
for  the  tube-and-disc  flameholder.  Thus  in  order  to  better  identify 
key  regions  of  inconplete  combustion,  flames  of  high  and  low  combustion 
efficiency  will  be  probed  for  selected  geometries  and  inlet  conditions . 
Both  time  averaged  species  concentration  (NO^,  CD,  HC,  0^,  CO^)  and 
teuperature  can  be  obtained  using  a single  point  probe  vdiich  is 
positioned  both  axially  and  radially  in  the  combustor.  This  probing 
technique  has  already  proven  successful  in  determining  regions  of 
NO^  and  CD  emissions  in  disc  stabilized  flames  (Tuttle  et  al.,  1975) 
and  in  turbojet  combustors  (Mellor,  1976).  Identification  of 
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appropriate  regions  of  combustion  inefficiency  greatly  assists  in  the 
choice  of  length  scales,  velocities  and  tenqieratures  relevant  to  the 
characteristic  time  model. 

To  further  characterize  the  flow  field,  we  are  also  planning 
the  design  and  construction  of  a laser  Doppler  velocimeter  (LDV) . The 
LDV  is  an  optical  device  \diich  makes  use  of  the  frequency  shift  in 
monochromatic  light  scattered  by  a moving  object  to  provide  a non- 
intrusive  means  of  measuring  local  flow  velocities.  Because  the 
frequency  shift  is  small  relative  to  the  frequency  of  the  incident 
light,  it  is  necessary  to  split  the  beam  and  use  one  part  as  a ref- 
erence. This  has  given  use  to  two  basic  LDV  designs:  the  reference 
beam  system,  in  which  only  one  of  the  beams  is  scattered  while  the 
other  is  focused  directly  onto  the  photomultiplier  tube;  and  the 
differential  LDV,  in  which  both  beams  are  focused  to  a point  in  the 
flow  and  scattered  light  from  both  beams  is  collected  and  processed. 
Because  the  output  of  the  reference  beam  LDV  is  dependent  on  the 
angle  at  which  the  scattered  light  is  collected  and  the  device  itself 
is  difficult  to  align,  we  will  make  use  of  the  differential  design. 

To  insure  the  presence  of  adequate  numbers  of  scattering  centers 
within  the  flow,  it  will  probably  be  necessary  to  introduce  artificial 
seeding.  In  so  doing,  however,  care  will  be  taken  to  insure  that  the 
particles  introduced  will  be  small  enough  to  follow  the  fluid 
accurately . 

Once  the  device  is  built  we  plan  to  make  detailed  maps  of 
the  combustor  flow  field  in  both  hot  and  cold  nonreacting  flows  and 
in  various  flames.  Though  the  number  of  axial  points  at  which  we 
make  measurements  in  the  test  section  is  limited  by  the  need  for 
optical  access  to  the  flow  to  the  five  window  positions,  the  radial 
resolution  of  the  device  should  allow  us  two-component  velocity  data 
at  twenty  or  more  points  at  each  such  axial  position.  By  using  our 
microcomputer  operated  data  acquisition  system  to  collect  and  reduce 
the  data,  we  will  be  able  to  take  enough  measurements  at  each  point 
in  a short  enougli  period  of  time  to  insure  a statistically  valid 
result. 
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These  detailed  velocity  magnitude  and  direction  measurements 
coipled  with  the  species  concentration  and  tenperature  profiles 
mentioned  earlier  should  aid  in  the  selection  of  scaling  parameters  for 
the  characteristic  time  model.  Furthermore,  these  detailed  measurements 
will  not  only  help  refine  our  model  but  will  also  provide  a data  base 
for  conparison  with  the  modular  modeling  techniques  propnosed  by 
Edelman  (1977). 
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Appendix  A 
Data  Listing 

All  of  the  combustion  efficiency  and  lean  blowoff  data 
obtained  for  the  simplified  prevaporizing/premixing  combustor  of 
Fig.  3-3  are  listed  in  this  section.  The  computer  print  out 
includes  nozzle  size,  fuel  type,  (starting  with  Jet  A),  geometry, 
pressure,  inlet  temperature,  air  flow  rate  and  overall  equivalence 
ratio.  (Note  that  the  equivalence  ratio  listed  here  is  an  apparent 
(j)  and  does  not  include  vitiation  effects) . In  addition  to  inlet 
conditions,  the  theta  parameter  of  Lefebvre  (1966)  has  been  computed 
for  each  data  point  as  well  as  the  initial  fuel  spray  SMD  and  the 
droplet  evaporation  time.  Parameters  in  are  evaluated  in  the 
same  manner  as  did  Tuttle  et  al.  (1976). 

Both  the  blowoff  and  combustion  efficiency  ejqjeriments  were 
performed  at  the  same  time  since  it  was  convenient  to  first  obtain 
stable  combustor  operation,  measure  the  emissions  and  then  slowly 
decrease  the  fuel  flow  rate  to  the  lean  limit.  The  comment  column 
indicates  whether  the  data  point  is  for  blowoff  or  combustion 
efficiency.  If  a blowoff  point  has  been  obtained,  all  emissions  and 
combustion  efficiency  values  are  printed  as  zeros.  However,  for  the 
combustion  efficiency  data,  a NO  El  value  of  zero  or  * indicates 
sanpling  problems  usually  associated  with  the  interaction  of  unbumed 
hydrocarbons  with  NO^  in  the  gas  sanT)ling  system  (Tuttle  et  al.,  1974). 
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Appendix  B 
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